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ABSTRACT 

He  discuss  and  interpret  a  theory  developed  by  Kreiss  and  others  for 
studying  the  suitability  of  boundary  conditions  for  linear  hyperbolic  systems 
of  partial  differential  equations.  The  existing  theory  is  extremely 
technical. 

The  present  discussion  is  based  on  the  characteristic  variety  of  the 
system.  The  concept  of  characteristic  variety  leads  to  : 

( 1 )  a  physical  interpretation  of  the  theory  in  terms  of  wave  propagation , 

and 

(2)  a  physical  and  geometrical  method  for  visualising  the  algebraic 
structure  of  the  system.  The  great  complexity  of  the  theory  is  caused  by 
certain  aspects  of  this  structure. 

We  also  point  out  connections  between  the  above  work  and  a  corresponding 
theory  regarding  the  stability  of  finite  difference  approximations. 
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SIGNIFICANCE  AND  EXPLANATION 


This  paper  is  concerned  with  boundary  conditions  for  linear  hyperbolic 
systems  of  partial  differential  equations  in  several  space  dimensions .  For 
various  reasons,  it  can  be  difficult  to  determine  whether  a  given  boundary 
condition  is  suitable  for  a  given  system.  The  existing  theory  which  deals 
with  this  issue  is  quite  complicated  and  algebraic.  In  the  present  paper  we 
describe  how  this  theory  can  be  interpreted  in  terms  of  concepts  associated 
with  wave  propagation.  We  also  mention  connections  with  an  analogous  theory 
regarding  the  stability  of  finite  difference  approximations. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 
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1 .  Introdnction. 

In  this  paper  we  disease  boundary  conditions  for  linear  hyperbolic 
partial  differential  equations  in  several  space  dimensions.  Our  goal  is  to 
provide  some  interpretations  and  explanations  of  a  theory  developed  by 
Kreiss  and  others  (e.g..  [10] , [14] , [151 , [20] , 121] )  for  determining  whether  a 
given  set  of  conditions  is  suitable  for  a  given  equation.  The  theory  is  of 
fundamental  importance.  However,  it  is  extremely  complicated,  and  its 
physical  interpretation  is  not  immediately  apparent.  In  the  present  section 
we  describe  the  context  of  the  theory  and  outline  the  interpretations  which 
we  will  give  regarding  it.  At  the  end  of  the  section  we  point  out  some 
relationships  between  this  theory  and  the  study  of  finite  difference 
approximations. 

Examples  of  hyperbolic  equations  include  the  Euler  equations  of  gas 
dynamics,  the  shallow  water  equations.  Maxwell's  equations,  equations  of 
magnetohydrodynamica,  and  the  classical  wave  equation.  Except  for  the  wave 
equation,  these  examples  are  systems  of  first-order  equations.  The  first 
three  esses  will  be  discussed  in  Section  4. 

The  theory  to  be  discussed  here  deals  with  linear  first-order  systems. 
Equations  of  higher  order  can  be  reduced  to  first-order  systems  by  standard 
techniques  (e.g.,  John  [9],  Taylor  [26]).  Of  the  above  examples.  Maxwell's 
equations  and  the  wave  equation  are  linear.  The  theory  is  applicable  to 
linearized  versions  of  the  others. 

The  unknown,  dependent  variables  in  the  problems  of  interest  are 
functions  of  space  and  time,  e.g.,  u  =  u(x,t)  where  x  8  Rn  .  The 
spatial  variable  x  is  typically  confined  to  a  subdomain  of  Ra  .  For 
example,  a  fluid  may  flow  in  a  region  which  is  bounded  by  a  solid  wall.  In 
other  problems  the  spatial  domain  of  interest  may  be  bounded,  at  least  in 
part,  by  an  open,  artificial  boundary.  These  boundaries  are  sometimes 
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introduced  in  order  to  limit  the  extent  of  a  problem  so  that  a  numerical 
computation  of  the  solution  can  be  made  feasible.  This  is  the  case  with 
limited-area  problems  in  oceanography  and  meteorology.  Another  such 
situation  is  the  modeling  of  a  fluid  flow  in  an  exterior  domain. 

In  order  to  determine  a  unique  solution  to  the  problem,  it  is  necessary 
to  specify  values  of  the  solution  at  some  initial  time,  and  it  is  generally 
necessary  also  to  impose  conditions  on  the  solution  at  the  boundary.  The 
problem  thus  becomes  an  initial-boundary  value  problem  (IBVP).  In  some 
cases  the  correct  boundary  conditions  can  be  found  easily  from  physical 
considerations.  At  a  solid  wall  which  bounds  the  flow  of  a  fluid,  for 
example,  one  sets  the  normal  component  of  the  fluid  velocity  equal  to  zero 
(if  effects  of  viscosity  are  to  be  considered,  the  tangential  component  must 
also  vanish).  In  other  situations  the  choice  of  boundary  conditions  is  not 
as  obvious.  This  is  the  case  with  artificial  boundaries,  which  do  not 
correspond  to  anything  physical. 

In  general,  it  is  necessary  to  be  careful  when  choosing  boundary 
conditions  for  a  hyperbolic  system.  This  can  be  seen  most  easily  in  the 
case  of  first-order  systems  in  one  space  dimension,  which  we  will  discuss  in 
Section  2.  There  we  will  show  that  various  portions  of  the  solution 
represent  traveling  waves.  It  will  be  apparent  that  any  acceptable  boundary 
condition  must  prescribe  the  behavior  of  the  waves  which  are  coming  into  the 
spatial  domain,  but  is  must  not  affect  the  waves  which  are  leaving. 

Examples  of  suitable  boundary  conditions  include  reflection  conditions, 
which  describe  incoming  waves  in  terms  of  outgoing  waves. 

In  several  space  dimensions  the  situation  is  more  complicated  than  in 
one  dimension.  In  this  case,  it  is  not  nearly  as  easy  to  identify 
"  incoming* '  and  ''outgoing'*  components  of  the  solution.  There  may  also  be 
waves  which  move  tangent  to  the  boundary,  and  it  may  not  be  clear  what  a 
boundary  condition  should  say  about  them.  These  difficulties  are  among  the 
principal  topics  to  be  discussed  in  the  present  paper. 

Because  of  the  above  problems,  it  may  be  difficult  to  determine  whether 
a  given  boundary  condition  is  suitable  for  a  given  equation  in  several 
dimensions.  One  approach  to  this  question  is  given  by  the  "energy  method" 
(e.g.,  Friedrichs  [4],  Courant  and  Hilbert,  vol.II  [3]).  This  method  gives 


criteria  which  are  sufficient  for  a  boundary  condition  to  yield  a  well-posed 
problem,  i.e.,  a  problem  which  admits  a  unique  solution  depending 
continuously  on  the  prescribed  data.  (See  Section  3.1  for  a  more  precise 
definition  of  ''well-posed".)  The  criteria  are  not  necessary;  if  the 
method  fails  to  show  that  a  boundary  condition  is  acceptable,  this  may  be 
dne  either  to  a  defect  in  the  condition  or  a  defect  in  the  method. 

An  alternate  approach  yielding  more  precise  information  ("normal  mode 
analysis")  has  therefore  been  developed.  The  theory  includes  work  by 
Kreiss  [10].  Sakamoto  [23].  Rauch  [21],  Ralston  [20],  Majda  and  Osher  [14], 
and  Michelson  [15]  .  The  theory  gives  criteria  which  are  essentially 
necessary  and  sufficient  for  a  boundary  condition  to  yield  a  well-posed 
problem. 

The  theory  is  extremely  complicated  and  algebraic;  its  physical 
meaning  tends  to  get  bnried  by  lengthy  and  detailed  studies  of  various 
matrices.  The  principal  purpose  of  the  present  paper  is  to  examine  this 
work  from  the  viewpoint  of  wave  propagation.  We  will  show  that  to  a  great 
extent  the  theory  has  physical  effects  that  one  would  expect  from  the 
discussion  of  incoming  and  outgoing  wavea  given  above.  The  concept  of  group 
velocity  plays  a  major  role  in  the  discussion.  The  great  complexity  of  the 
theory  is  due  mainly  to  the  Fourier  components  of  the  solution  which 
correspond  to  waves  traveling  tangent  to  the  boundary. 

We  now  outline  the  contents  of  this  paper.  In  Section  2  we  review  the 
standard  treatment  of  systems  in  one  space  dimension.  In  Section  3  we 
describe  the  motivation  which  is  usually  given  for  the  multi-dimensional 
theory.  Section  3  concludes  with  a  literature  survey. 

The  interpretations  to  be  given  here  make  extensive  use  of  the 
structure  of  the  characteristic  variety  of  the  system,  i.e.,  the 
high-frequency  part  of  the  "dispersion"  relation.  In  Section  4  we 
describe  this  structure  for  some  examples,  and  we  point  out  that  these 
special  cases  model  crucial  features  of  more  general  systems.  For  the  sake 
of  definiteness  and  clarity,  the  discussion  in  the  remaining  sections  is 
given  in  terms  of  the  examples.  In  Proposition  4.2  (Section  4.4)  we 
describe  a  process  by  which  incoming  and  outgoing  modes  are  labelled 
algebraically  in  the  theory  of  well-posedness. 
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In  Section  S  we  give  sone  interpretations  of  the  ' 'uni for*  Kreiss 
condition* *  (D.K.C.),  the  slgebrsic  criterion  which  chsrseterixes  admissible 
boundary  conditions.  This  criterion  is  nsed  to  solve  for  ’'incoming*' 
components  in  terms  of  "outgoing"  components  and  boundary  data.  In  this 
section,  we  also  consider  some  weaker  forms  of  well-posedness  which  are 
encountered  when  the  U.K.C.  is  not  quite  satisfied.  In  these  cases  there 
are  certain  uniformities  which  do  not  hold  as  the  frequencies  of  the  Fourier 
components  of  the  solution  tend  to  infinity,  so  that  certain  reflection 
coefficients  may  tend  to  infinity. 

In  Section  6  we  use  the  framework  of  Sections  4  and  5  to  discuss  some 
of  the  main  ideas  in  the  proofs  by  Kreiss  and  others  which  show 
well-posedness  of  the  initial-boundary  value  problem. 

Sections  4,5,  and  6  form  the  core  of  this  paper.  The  reader  who  is  in 
a  hurry  may  find  it  worthwhile  to  skim  lightly  over  Section  2  and  Section 
3.2. 

The  theory  discussed  here  is  closely  related  to  a  theory  which  deals 
with  the  stability  of  finite  difference  approximations  to  initial-boundary 
value  problems.  The  stability  theory  includes  work  by  Osher  [18], [19], 
Gustafsson,  Kreiss.  and  Sundstrom  [5],  and  Nicholson  [16].  Trefethen 
[27], [29]  has  recently  studied  this  work  from  the  viewpoint  of  wave 
propagation  and  has  reached  some  conclusions  which  are  analogous  to  some 
ideas  expressed  in  the  present  paper.  From  time  to  time  we  will  point  out 
some  similarities  between  the  well-posedness  theory  and  the  stability 
theory. 
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2.  Systems  in  one  space  dimension 


We  now  rewiew  the  situation  for  systems  in  one  space  dimension.  In 
this  ease  the  problem  of  finding  suitable  boundary  conditions  is  fairly 
straightforward. 

We  consider  the  system 


(2.1)  ut  -  Aux  +  Cu  , 


where  u(x,t)  is  a  vector  having  n  components,  and  A(x,t)  and  C(x.t) 
are  n  x  n  matrices.  For  simplicity,  we  assume  that  (2.1)  is  defined  for 
x  >  0  and  t  >  0.  If  a  system  is  defined  on  a  spatial  domain  with  two 
boundaries,  e.g.,  a  <  x  <  b  ,  then  each  boundary  can  be  treated  separately 
in  the  manner  to  be  described  below. 

The  system  (2.1)  is  assumed  to  be  hyperbolic,  i.e.,  A  has  real 
eigenvalues  and  a  complete  set  of  real  eigenvectors.  This  assumption 
enables  one  to  simplify  the  form  of  the  system.  It  implies  that  there 
exists  a  nonsingular  matrix  Q  such  that  Q-1AQ  *  diag  ^....Xg). 

The  system  (2,1)  can  then  be  written  as 


(Q-1u)t  -  (Q-1AQ) (Q-1u) x  +  (Q-1C  +  Q’J  -  Q_1AQQ_J)u  , 


or 


(2.2) 


vx  +  Dv 


where  v  =  Q~*u  . 

For  the  sake  of  simplicity,  we  first  consider  the  case  D  =  0  .  Under 

this  assumption  each  equation  in  (2.2)  has  the  form 

dv.  dv, 

— 1  -  X,  — 1  »  0 
dt  ■*  dx 

and  is  an  ordinary  differential  equation  for  Vj  along  the 
characteristic  curves  defined  by  dx/dt  ■  -Xj .  (See  Figure  2.1)  The 
components  Vj  are  constant  along  the  corresponding  characteristic 
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curves  and  can  thus  be  regarded  as  traveling  waves  which  move  at  the 
characteristic  velocities  dz/dt  “  -Xj . 


A.  <  0  A  >  0 

J  j 


Figure  2.1.  Characteristics  for  (2.2). 

Initial  values  for  the  ordinary  differential  equations  for  Vj  are 
provided  by  the  values  of  Vj  at  t  **  0  and  also,  when  Xj  <  0  , 
by  values  of  Vj  at  the  boundary  z  =  0  .  It  is  therefore  necessary  to 
prescribe  an  initial  condition 

v(z,0)  *  f(z)  ,  z  >  0  , 

where  f  is  a  given  function,  together  with  a  boundary  condition  which 
defines  values  for  the  vj  corresponding  to  negative  Xj  ,  i.e., 
the  incoming  components  of  v  .  It  is  not  permissible  to  prescribe  values 
for  outgoing  components  at  z  ■  0;  otherwise,  the  boundary  conditions 
could  contradict  the  effect  of  the  initial  condition  and  thereby  make  it 
impossible  for  a  solution  to  ezist.  The  boundary  conditions  must  therefore 
fit  the  general  form 

(2.3)  vI(0, t)  =  SvII(0, t)  +  g(t>  ,  t  >  0  , 

where  v*  and  v**  are  vectors  consisting  of  the  incoming  and 
outgoing  components,  respectively.  Here  g  is  a  given  function,  and  S  is 
a  rectangular  matriz  of  appropriate  dimensions  which  governs  reflections  at 
the  boundary. 

When  D  #  0  ,  the  equations  in  (2.2)  are  coupled  together.  In  this 
case  the  ezistence  of  the  solution  can  be  shown  via  an  iteration  of  Picard 
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type.  See,  e.g.,  Courant  end  Hilbert  t3].  The  boundary  conditione  aust 
satisfy  the  saae  criteria  as  before;  it  is  necessary  to  prescribe  values 
for  incoaing  coaponents,  and  it  is  not  peraissible  to  prescribe  values  for 
outgoing  coaponents.  The  identification  of  the  incoaing  and  outgoing 
coaponents  depends  only  on  the  leading  order  part  of  the  systea,  i.e.,  is 
independent  of  D  . 

It  is  often  of  interest  to  know  that  the  solution  depends  continuously 
on  the  prescribed  data  f  and  g.  In  this  case  the  continuous  dependence 
follows  froa  the  fact  that  solutions  to  ordinary  differential  equations  are 
continuous  functions  of  their  initial  data.  (Also  see  (3.5)  and  the 
associated  discussion) 
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3.  Maitaiiaa  isx  ihs. 


theory 


ft  now  begin  the  disenssion  of  problem  in  several  space  dimensions. 

In  this  case  it  is  generally  not  possible  to  construct  and  analyze  the 
solution  by  using  characteristic  curves.  As  noted  in  the  Introduction, 
additional  problems  may  be  caused  by  the  possibility  of  waves  moving  tangent 
to  the  boundary.  One  can  therefore  expect  this  case  to  be  more  complicated 
than  the  previous  one. 

In  the  present  section  we  describe  the  motivation  which  is  usually 
given  for  the  approach  developed  by  Kreiss  and  others  for  studying  the 
multi-dimensional  case.  In  Section  3.1  we  make  some  preliminary  remarks, 
and  in  Section  3.2  we  derive  a  condition  which  is  necessary  for  the  problem 
to  be  well-posed.  In  Section  3.3  we  state  a  stronger  condition  which  can  be 
shown  to  be  sufficient.  The  effects  of  the  stronger  assumptions  will  be 
discussed  in  later  sections,  mainly  Section  5.  The  proof  of  sufficiency  is 
the  main  source  of  difficulty  in  the  theory  and  bears  little  resemblance  to 
the  derivation  of  the  necessary  condition.  In  Section  3.4  we  survey  some  of 
the  literature  which  deals  with  the  subject. 

3.1  Preliminaries 

A  fairly  general  form  for  linear  hyperbolic  systems  is 

m 

(3.1)  Uj.  =  Au,  +  7  B-  —  +  Cu  +  F 

‘  *  J-l  J*7j 

(S»e  the  examples  in  Section  4.)  We  consider  this  system  for  t  >  0  ,  and 
for  reasons  given  below  we  assume  that  the  spatial  domain  is  defined  by 
x  >  0  and  y  e  R®,  where  ■  2  1<  In  (3.1)  u(x,y,t)  and  F(x,y,t) 
are  vectors  having  n  components,  and  A,Bj,  and  C  are  n  x  n 
matrices  which  we  assume  to  be  smooth  functions  of  x,y,  and  t. 

In  the  theory  which  has  been  developed,  the  system  is  assumed  to  be 
either  strictly  hyperbolic  or  syauetric  hyperbolic.  In  the  former  case  the 
awtrices  oA  +  £  UjBj  have  real  distinct  eigenvalues  for  all  real  a, 

-j.  ...  «B  for  which  lor  I  +  lul  /  0  ,  In  the  latter  case  A  and  the 


Bj  are  hermitian.  (See  Section  4.1  regarding  ' ' symmetrizable ' *  systems.) 

The  spatial  domain  [(>,•)  z  R®  has  been  chosen  for  the  sake  of 
simplicity.  If  one  is  considering  a  system  defined  on  a  spatial  region 
which  does  not  have  this  form  but  still  has  a  smooth  boundary,  then  one  can 
localize  the  problem  with  a  partition  of  nnity  and  then  map  each  boundary 
portion  into  the  boundary  of  the  half-space  [0,*)  x  R®  .  In  the  new 
coordinates  the  problem  will  have  the  form  given  above. 

Since  (3.1)  is  hyperbolic,  the  matrix  A  can  be  assumed  to  be  diagonal 
with  real  eigenvalues;  otherwise  one  can  find  a  similarity  transformation 
which  makes  it  diagonal  and  then  adopt  the  corresponding  change  of  dependent 
variable.  In  much  of  the  theory  A  is  also  assumed  nonsingular,  i.e.,  the 
boundary  z  =  0  is  noncharacteristic,  (cf.  Section  4.3)  Unless  otherwise 
stated,  we  will  assume  that  this  is  the  case  and  that  the  elements  of  A 
are  arranged  so  that  A  has  the  form 


where  A*  <  0  and  A*1  >  0  . 

In  analogy  with  the  one-dimensional  case,  we  prescribe  an  initial 
condition 

(3.3)  u(x,y,0)  =  f (x,y)  , 

where  f  is  a  given  function,  together  with  a  boundary  condition  of  the 
form 

(3.4)  u^(0,y, t)  =  S  uII(0,y,t)  +  g(y,t)  . 

Here  g  is  given  function,  S  is  a  rectangular  matrix,  and 
n  *  ((®^)^,  (u11)*)*  .  The  components  of  u^  and  u^  correspond  to 
the  blocks  A*  and  A*^  ,  respectively,  in  (3.2).  Some  remarks 
about  the  form  of  (3.4)  are  given  at  the  end  of  Section  5.1. 

The  problem  at  hand  is  to  determine  whether  the  boundary  condition 

(3.4)  is  appropriate  for  the  system  (3.1).  In  the  one-dimensional  case  the 


answer  is  obvious,  but  in  the  present  case  it  nay  be  necessary  to  place  some 
restrictions  on  the  matriz  S  .  The  effect  of  the  theory  being  described 
here  is  to  identify  restrictions  which  are  necessary  and  sufficient  for  the 
initial-boundary  value  problem  to  be  well-posed. 

By  ''well-posed''  we  mean  that  for  arbitrary  F.f.g  in  suitable 
function  spaces,  the  problesi  (3.1),  (3.3),  (3.4)  admits  a  unique  solution, 
and,  furthermore,  it  is  possible  to  estimate  the  solution  in  terms  of 
F, f,g.  The  latter  condition  is  equivalent  to  the  continuous  dependence  of 
the  solution  on  the  prescribed  data.  A  typical  estimate  ("energy 
estimate")  has  the  form 


(3.5)  ll*l,9*to.«J  *  ll»ll3(!xl..tl  * 

where  K  is  independent  of  u,f,F,  and  g  (e.g.,  Majda  and  Osher  [14], 

Rauch  [21]).  Q  denotes  the  spatial  domain  z  >  0  ,  and  the  norms  are 
2 

weighted  L  norms  or  Sobolev  norms  on  the  regions  indicated  by  the 
subscripts.  See  also  the  estimate  (5.25). 


3.2.  A  necessary  condition  for  wcll-nosedness. 

Ve  now  develop  the  necessary  condition  which  was  mentioned  in  the 
introduction  to  Section  3.  This  condition  is  based  on  certain  families  of 
ezponential  solutions  of  the  differential  equation  which  cannot  possibly 
satisfy  the  energy  estimate  (3.5),  or  any  similar  estimate.  If  these 
functions  also  satisfy  the  boundary  conditions,  then  they  are  solutions  of 
the  initial-boundary  value  problem  (IBVP),  and  the  problem  must  be 
ill-posed.  It  is  therefore  necessary  to  guarantee  that  the  boundary 
conditions  ezclude  these  functions. 

The  necessary  condition  is  similar  to  the  Godunov-Ryabenkii  criterion 
for  the  stability  of  finite  difference  approzimations.  (See,  e.g.,  [5], [22], 
or  [28]). 

In  order  to  obtain  the  special  solutions,  we  assume  that  C  «■  0  ,  A, 

Bj  and  S  are  constant,  and  F  «  g  ■  0  .  The  IBVP  (3.1),  (3.3),  (3.4) 
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is  then 


(a)  u*  “  An 


_ 

‘  *  1  ^ 


,x>0,  y*R*»t>0 


(3.6) 


(b)  u(x,y,0)  ■  f(x.y) 


(c) 


Su 


II 


The  solntions  of  interest  ("normal  nodes' ')  are  elementary  waves 
having  the  form 


(3.7) 


n  -  4<x)  e1-*  +  st  . 


where  «  e  R®  and  s  is  complex  with  Be  s  >  0  Modes  of  this 
type  are  associated  with  Laplace  and  Fonrier  transforms  (see  Section  4),  but 
we  will  not  need  to  nse  these  transforms  at  the  present  time.  In  fact, 
their  nse  would  introduce  unnecessary  complications. 

If  det  A  ^  0,  (3.6) (a)  can  be  written  as 


(3.8)  ux  =  A-1(nt  -  [  B,  —  ) 

J  *7j 

When  (3.7)  is  substituted  into  (3.8)  ,  the  result  is  an  ordinary 
differential  equation  for  the  amplitude  function: 

m 

(3.9)  4*(x)  -  A-* (si  -  T  iw.B . )  4 

j-i  J  J 

-  m(«,s)4 

The  solntions  of  (3.9)  will  be  discussed  extensively  later.  For  the  moment 
it  suffices  to  say  that  there  are  certain  solntions  which  decay 
exponentially  as  x  •*  +  •  and  others  which  grow  exponentially.  Ve  also  note 
that  for  any  solution  n  of  the  form  (3.7),  the  related  functions 

(3.10)  ua(x.y.t)  -  4(ax)  eiaw*7  +  0,4 

are  also  solntions  of  (3.6) (a),  for  all  real  o. 
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A  necessary  condition  for  well-posedness  is  iaplied  by  the  following 
Proposition: 

Proposition  3.1.  Snppose  that  for  some  s  with  Se  s  >  0  and  sosie 
m  s  SF  ,  there  is  a  function  u  of  the  fora  (3.7)  which  satisfies 
the  following: 

(a)  4*0.  and  4  decays  exponentially  as  x  -»■  +•  ;  and 

(b)  n  satisfies  the  bonndary  condition  (3.6)(c),  i.e., 

(3.11)  u*  “  Su**  ,  for  x  «  0  . 

The  IBVP  (3.1),  (3.3),  (3.4)  aust  then  be  ill-posed. 

Proof.  (3.11)  iaplies  that  the  corresponding  na  also 
satisfy  the  bonndary  condition,  since  (3.11)  is  eqniwalent  to  4*(0)  ■ 
S4II(0).  The  nQ  are  thus  solutions  of  the  IBVP  (3.6),  which 
is  a  special  case  of  (3.1),  (3.3),  (3.4).  Ve  also  note  that  the  nQ 
have  finite  nora  with  respect  to  x  . 

We  now  show  that  the  faaily  (ua:  a  >  0}  rlol*te$  the  energy 
estinate  (3.5).  As  a  -*■  +-  ,  the  solutions  uQ  grow  at 

arbitrarily  high  exponential  rates  in  t  ,  since  Re  s  >  0  .  Equivalently, 
the  left-hand  ride  of  (3.5)  is  an  exponentially  increasing  function  of  o. 
This  is  not  the  case  for  the  right-hand  side.  In  this  case  F  and  g  are 
zero,  and  the  initial  values  are 

fa(x.y)  -  4 (ox)  eiaw*y  . 

Sobolev  noras  with  respect  to  x  and  y  introduce  derivatives  which  cause 
polynoaial  growth  in  a  ,  but  this  is  the  most  that  can  happen  on  the 
right-hand  side.  By  taking  a  sufficiently  large,  we  can  therefore 
conclude  that  no  constant  K  in  (3.5)  is  adequate.  The  definition  of 
well-posedness  given  earlier  is  thus  violated. 

Hersh  [6]  has  pointed  out  that  the  failure  of  the  energy  estiaate 
iaplies  that  existence  and/or  uniqueness  aust  also  fail  to  hold.  Bis 


12 


argument  is  based  oa  the  closed  graph  theorem. 

A  crucial  part  of  the  above  argosieat  is  the  fact  that  F  and  g  are 
zero.  If  F  and  g  are  large  whenever  the  solntion  is  large,  then  there 
is  no  problea.  Bat  in  the  case  considered  here,  it  is  possible  for  large 
eolations  to  appear  withont  sufficient  provocation  from  the  prescribed  data. 

One  any  object  that  the  nQ  do  not  have  finite  L  norms  or 
Sobolev  norms  with  respect  to  y  and  therefore  wonld  not  have  to  satisfy 
(3.5).  This  problem  can  be  remedied  by  truncating  the  nQ  .  Let  1 p 
be  a  smooth  function  of  y  which  is  equal  to  1  for  lyl  £.  1  and  is  equal 
to  zero  for  lyl  2  2  .  and  consider  the  functions 

u0  jt(x,y,t)  =  =  ^|/(^)^(ax)eia<,,' y  +  ast 

for  X  >  0  .  These  functions  satisfy  an  inhomogeneous  differential 
equation  and  hoaogeneous  boundary  conditions.  The  failure  of  the  energy 
estimate  is  shown  by  letting  a +«  .  We  omit  the  details. 


The  ideas  in  the  above  proof,  aside  from  the  truncation  in  y  ,  have 
been  credited  to  Agm on.  (cf.  Kreiss  110),  Agmon  II]) 

The  function  4  of  Proposition  3.1  is  often  said  to  be  the  solution 
of  an  eigenvalue  problem  with  eigenvalue  s,  since  it  satisfies 

d4 

si  =  A —  +  (  Ji«.B,)i 
dx  J  J 

(3.12)  ir(0)  -  SiI]C(0) 

4  e  L2(0.«) 


The  third  condition  may  be  regarded  as  a  boundary  condition  at  infinity. 
Proposition  3.1  can  then  be  stated  in  the  following  manner: 


Pronosit ion  3.1*  .  The  IBVP  (3.1),  (3.3),  (3.4)  is  ill- posed  if,  for 
some  <■>  e  R®,  the  problem  (3.12)  has  an  eigenvalue  s  with  Re  s  >  0. 


We  will  later  formulate  this  Proposition  as  an  algebraic  criterion 
which  is  necessary  for  well-posedness,  but  we  must  first  characterize  the 
eolations  of  the  system  (3.9). 
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Tbe  basic  theory  of  ordinary  differential  equations  implies  that  (3.9) 
has  n  linearly  independent  eolations  which  span  the 

set  of  all  solutions  of  (3.9).  The  can  be  taken  to  have  the  fora 

(3 .13)  (a)  4(x)  ■  e<x  w 

where  <  is  an  eigenvalue  of  M(w,s)  and  w  is  a  corresponding 
eigenvector,  or 

(3.13)  (b)  4(x)  »  e<xP(x)  . 

where  P  is  a  polynoaial  having  vector  coefficients.  The  latter  fora 
arises  when  <  is  associated  with  a  nontrivial  Jordan  block.  The  degree  of 
P  is  less  than  the  algebraic  anltiplicity  of  K  ,  and  P(0)  is  a 
generalized  eigenvector.  Ve  will  not  denote  explicitly  the  dependence  of 
4.  w,  and  P  on  s  and  «  . 

The  4j  give  rise  to  exponential  solntions  of  the  partial  differential 
equation  (3.6) (a)  via  the  relation  (3.7).  a  *  4(x)ei“*y  +  xt.  These  aodes 
will  be  denoted  by 

Ojix.y.t)  -  ♦j(*)ei-y  +  “ 

(3.14)  f  e<x  +  im'7  *  •*  w 

■  4 

or  e<x  +  ltt'r  +  #t  P(x) 

Ve  are  interested  in  solutions  of  this  type  which  have  finite  L2  noras 
or  Sobolev  noras  with  respect  to  x  on  the  interval  0  <  x  <  •  . 

Inforaation  about  these  aodes  is  given  by  the  following  Proposition  due  to 
Hersh  (6).  Analogous  properties  of  finite  difference  equations  are  given  in 
Leaaas  5.1  and  5.2  of  Gustafsson,  Ireiss,  and  Sundstroa  [5]. 

Proposition  3.2.  Let  l  denote  the  nuaber  of  negative  eigenvalues  of 
A.  If  te  s  >  0.  then  the  aatrix  M(v.s)  »  A~2(sl  -  l  iwjBj) 


(see  (3.9))  has  £  eigenvalues  <  with  negative  real  part  and  n  -  & 
eigenvalue*  with  positive  real  part. 

Proof.  Ve  first  note  that  there  are  no  purely  ina ginary  eigenvalues 
when  Re  s  >  0;  if  <  *  io  were  an  eigenvalue,  then 

iow  *  A~*(sl  -  l  iUjBj)w  , 

for  some  vector  w  j*  0,  and  therefore 

(ioA  +  £iWjBj)w  »  sw  . 

But  s  would  then  have  to  be  purely  imaginary,  since  the  system  is 
hyperbolic . 

Let  s  *  n  +  i£.  Ve  are  interested  in  the  signs  of  the  real  parts  of 
the  eigenvalues  of  N(w, s)  on  the  domain  n  >  0,  5  e  R,  »  e  R®. 

The  eigenvalues  are  continuous  functions  of  n,5,  and  u,  so  the 
real  part  of  each  eigenvalue  must  have  a  constant  sign  on  the  domain; 
otherwise  there  would  be  sets  on  which  Re  K  *  0  for  some  k  .  Ve  can 
therefore  count  the  signs  at  any  point  which  may  be  convenient.  If  n  *  1. 
5  =  0,  and  u  =  0,  then  H  »  A-1.  In  this  case  the  claimed  counting 
is  correct,  and  the  Proposition  follows. 

■The  Proposition  implies  that  it  of  the  functions  in  (3.13) 

have  finite  norm  on  the  interval  0  <  x  <  «  .  The  same  is  true  of  the 
solutions  Uj  in  (3.14),  for  fixed  y  and  t.  Arrange  indices  so  that 
these  functions  are  ^....^  and  uj,...  u^. 

Ve  can  now  formulate  an  algebraic  condition  which  is  necessary  for 
well-posedness.  According  to  Proposition  3.1,  it  is  necessary  to  prevent 
certain  solutions  u  of  (3.6) (a)  from  satisfying  the  homogeneous  boundary 
condition  (3.6) (c).  The  solutions  of  concern  are  linear  combinations  of 
Uj,...  u£*  an<*  the  corresponding  d's  are  linear  combinations  of 
.... 

The  boundary  condition  (3.6) (c)  can  be  written  as 
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»  0,  for  x  ■  0 
For  the  solutions  of  interest  this  condition  has  the  fora 
(3.15)  [I,  -SJIuj,...  u^l  c  -  0  .  for  x  -  0. 

where  c  is  sn  £  -  vector  whose  components  are  the  coefficients  in  the 
linear  coaibination.  The  representation  Uj  ■  4j(x)eiw  7  +  (gee  (3.14)) 
iaiplies  that  (3. IS)  is  equivalent  to 

(3.14)  [I.  -S] [4^(0) , . . .  4£(0)]  c  ■  0 

The  linear  independence  of  4j,...  4^  implies  that  a  linear 
combination  of  u^,...  n^  is  xero  if  and  only  if  c  “  0.  Ve  therefore 
want  (3.16)  to  have  no  nonzero  solutions  c,  which  means  that  the  matrix 

(3.17)  N(«.s)  -  U.-S)[41(0)....4fc(0)] 
must  be  nonsingular. 

Ve  should  note  that  N(w, s)  is  square;  the  number  of  rows  in  the 
boundary  condition  (3.6) (c)  is  equal  to  the  number  of  negative  eigenvalues 
of  A,  which  is  i  .  (See  (3. 2), (3. 4),  and  Proposition  3.2.)  Ve  also  note 
that  N  really  does  depend  on  «  and  s;  according  to  (3.13)  and  the 
associated  explanation,  the  vectors  4^(0) ,.. .4^(0)  are  eigenvectors 
and  generalised  eigenvectors  of  N(«, s).  They  correspond  to  the 
eigenvalues  which  have  negative  real  part.  For  an  example  of  N(«, s),  see 
(3.23). 

From  the  above  discussion  we  can  conclude  the  following: 

Proposition  3.3.  A  necessary  condition  for  well-posedness  is  that 

(3.18)  det  N(m.s)  4  0 
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for  oil  •  i  I*  and  all  complex  a  with  Be  a  >  0. 


3.3.  Xfcq  toff  1c lent  goaflitlQB 

We  nov  introduce  a  stronger  version  of  (3.18)  wbich  can  be  shown  to  be 
anfficient  for  well-posedness.  The  stronger  condition  is  a  an i fora  veraion 
of  (3.18),  and  in  at  least  many  cases,  it  amounts  to  a  requirement  that 

(3.18)  hold  for  Re  a  2  0,  rather  than  jnat  for  Re  a  >  0.  The  effecta  of 
the  atronger  assumptions  will  be  discussed  in  later  sections. 

In  order  to  state  the  sufficient  conditon,  it  will  be  necessary  to 
introduce  some  normalizations.  For  any  (u,s)  with  Re  a  >  0,  choose  a 
basis  for  the  &  -  dimensional  vector  apace  of  decaying  solutions  of  the 
system  4X  ■  M(w,s)4  ( i.e. . (3 .9) ) ,  and  from  this  basis  produce  a 
new  basis  which  is  orthonormal  at  x  =  0.  For  example,  one  could  choose  the 
particular  functions  4^»...  4^  discussed  earlier  (see  (3.13)  and 
the -discussion  after  Proposition  3.2)  and  then  perform  the  Gram  -  Schmidt 
orthogonal izat ion  process  on  the  initial  values  4^(0),...  4^(0). 

We  have  noted  that  these  initial  values  are  eigenvectors  and  generalized 
eigenvectors  of  M(w,s).  For  reasons  given  below,  we  will  not  limit 
attention  to  the  particular  basis  4j,... 4^  . 

Denote  the  orthonormal  initial  values  by  *)»••  (»,s) ; 

and  let 

(3.19)  N(m.s)  -  [I,-Smi«,s)...4>jl(«,s)] 

(cf.(3.17)).  The  sufficient  condition  is  the  following: 

Uniform  Erelss  Condition  (U.t.C.) .  There  exists  b  >  0  such  that 

(3.20)  Idet  N(u.s) I  2  6 

for  all  »  e  H*  and  all  complex  s  with  Re  s  >  0. 
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This  i*  essentially  the  formulation  used  by  Majda  and  Other  [14]  end 
Michelson  [IS].  An  alternate  formulation  will  be  discussed  later.  Ve  make 
several  comments  regarding  the  form  of  (3.20): 

(1)  The  earlier,  necessary  condition  (3.18)  did  not  involve  any 
normalization  of  the  vectors  ^(0),...  4j^0),  but  some  sort  of 
normalization  is  necessary  if  the  uniform  condition  (3.20)  is  to  hold. 
Otherwise,  we  coold  replace  ijij  with  etjjj,  with  e  small. 

(2)  The  U.I.C.  does  not  depend  on  the  choice  of  orthonormal  basis 

I|i1(w,s),...^(*>,s).  If  I is  another  such  basis,  then 

[^....^1  -  [i^,...^]Q. 

where  0  is  an  l  x  l  matrix.  It  is  easy  to  check  that  Q  is  unitary. 
Then, 

ldet([I,-smi....5jJ])l  -  det([I,-Sl['l>1...>£]Q)l 

-  Idet(NQ) I  =  Idet  Nl  *  Idet  Ql 

-  Idet  Nl 

(3)  It  would  suffice  to  impose  (3.20)  merely  for  Iwl^  +  IsP  *  1 
(and  Re  s  >  0);  in  each  direction  in  the  (m,s)  space  it  is  possible 

to  choose  vectors  ^j(u,s)  which  are  homogeneous  of  degree  zero  in 

(**, s).  For  example,  this  can  be  accomplished  by  observing  that 

M(«, s)  “  A”*(sl  -  ^ifeijBj)  is  homogeneous  of  degree  one  and  that  solutions 

of  «  W)  can  be  written  in  terms  of  solutions  for  l*»P  +  Isl^  *  1 

by  scaling  x  appropriately. 

(4)  The  FKC  applies  to  systems  which  have  variable  coefficients.  For 
such  systems  the  matrix  M(w,s)  is  also  a  function  of  (x,y,t).  In  this 
case  the  above  analysis  is  applied  to  the  constant-coefficient  problems 
obtained  by  freezing  coefficients  at  boundary  points  (0,y,t),  and  the  UKC  is 
then  required  to  hold  uniformly  in  (y.t).  The  calculus  of 
pseudo-differential  operators  is  used  extensively  in  the  treatment  of 
variable-coefficient  problems.  (See  Section  6). 

Ve  now  discuss  an  alternate  formulation  of  the  U.K.C.  The  earlier 
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condition  (3.18) ,  dot  N(w, s)  t  0  for  Re  s  >  0,  is  equivalent  to 
requiring  det  N(w, s)  ft  0  for  Re  s  >  0.  This  follows  easily  froa 
s  change-of-bssis  argnnent  which  reseables  that  of  comnent  (2)  above.  Ve 
now  show  that  in  aany  cases  the  U.K.C.  is  equivalent  to  requiring 
det  N(w, s)  f1  0  for  Re  s  2  0,  not  just  for  Re  s  >  0.  The  U.K.C.  is  often 
described  in  this  aanner  (e.g.,  in  13] ,110]). 

The  equivalence  will  follow  froa  soae  continuity  arguaents.  In  at 
least  aany  cases  (see  Section  6),  it  is  possible  to  choose  the  basis  vectors 
^(u.  s) , . .  (w,  s)  so  that  they  are  piecewise  continuous  functions  of  (a.s). 

The  continuity  ashes  it  possible  to  extend  the  definition  of  N(u,s)  to 
Re  s  «  0.  Ve  note  that  the  description  of  N(w, s)  given  earlier  is  not  valid 
in  this  case;  when  Re  s  ■  0  the  systea  ix  **  Mi  does  not,  in  general, 
have  £  linearly  independent  solutions  which  tend  to  zero  as  x  -*  +  • 

(cf.  Proposition  3.2  and  Section  4.3).  The  equivalence  aentioned  above  now 
follows  froa  the  piecewise  continuity  of  N(w, s)  and  the  coapactness 
of  the  set  ((w,s):  l«l2  +  Isl2  •  1,  Re  s  2  0) . 

3.4  A  survey  of  soae  of  the  literature  on  the  sublect. 

Hersh  [6]  studied  first-order  hyperbolic  systeas  which  have  constant 
coefficients  and  which  are  defined  on  a  half-space.  He  showed  existence  and 
uniqueness  of  solutions  subject  to  the  condition  (3.18).  Solutions  were 
constructed  via  Fourier  transforas  and  Laplace  transforas. 

Kreiss  [10]  used  an  approach  which  is  applicable  to 
variable-coefficient  probleas  defined  on  arbitrary  spatial  domains  having 
saooth  boundaries.  Here  the  aain  idea  is  to  construct  s  "  syaaetrizer' ' ,  a 
pseudo-differential  operator  which  enables  one  to  obtain  the  necessary 
energy  estiaate.  This  is  the  approach  which  will  be  discussed  later  in  this 
paper.  The  construction  is  based  on  the  condition  (3.20).  Kreiss  assumed 
that  the  system  is  strictly  hyperbolic  and  that  the  matrix  A  in  (3.1)  is 
nonsingular,  i.e.,  the  boundary  x  ■  0  is  noncharacteristic.  He  also 
assuaed  that  the  initial  data  are  zero. 

Sakaaoto  [23]  obtained  analogous  results  for  scalar  hyperbolic 
equations  of  higher  order. 


Ralston  [20]  provided  an  alternate  method  for  treating  soae  technical 
points  in  Kreiss'  construction.  Bis  work  extends  Kreiss’  resnlts  to  systems 
having  coaplex  coefficients. 

Rauch  [21]  considered  the  case  of  nonxero  initial  data  and  proved 
regularity  estimates  for  the  solution. 

Majda  and  Osher  [14]  generalised  the  theory  to  the  case  where  the 
boundary  is  "uniformly  characteristic".  This  aeans  that  the  aatrix  A  ■ 
A(x,y,t)  in  (3.1)  has  constant  rank  <  n  for  all  (x.y.t)  in  a  neighborhood 
of  the  boundary  x  *  0  .  They  assumed  that  the  system  is  syaaetrixable 
hyperbolic. 

Michelson  [IS]  used  a  theory  of  analytic  matrices  to  simplify  some 
aspects  of  the  uniformly  characteristic  case. 

Strikwerda  [25]  studied  "incompletely  parabolic"  systems  from  a  point 
of  view  which  is  similar  to  that  of  Kreiss.  et.al. 

Sarason  and  Smoller  [24]  used  methods  of  geometrical  optics  to  study 
the  behavior  of  hyperbolio  systems  in  spatial  regions  with  corners.  They 
considered  the  case  of  rays  which  are  not  tangent  to  the  boundary. 

Majda  [13]  used  Kreiss  symmetrizers  to  study  the  linearized  stability 
of  multi-dimensional  shock  fronts. 

The  well-posedness  theory  discussed  here  is  closely  related  to  a  theory 
for  the  stability  of  finite  difference  approximations  to  initial-boundary 
value  problems  for  hyperbolic  systems.  Gustafsson,  Kreiss.  and  Sundstrom 
[5]  showed  that  a  condition  analogous  to  the  U.K.C.  implies  the  stability  of 
difference  approximations  in  one  space  dimension.  Osher  [18].  [19]  had 
earlier  used  different  techniques  to  obtain  stability  results  for  a 
restricted  class  of  methods.  Michelson  [16]  has  obtained  results  for 
dissipative  approximations  in  several  dimensions  when  the  boundary  is 
noncharacteristic  and  the  initial  data  are  zero.  A  more  extensive  survey  of 
the  literature  is  given  by  Trefethen  [27]. 
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4 


IlXla&U  cxaaolcs. 


la  this  section  we  disease  various  Betters  related  to  the 
''characteristic  variety"  of  the  systea  (3.1).  The  discussion  will  give 
useful  inforastion  about  the  nature  of  the  well-posedness  theory  and  its 
physical  interpretation. 

Before  we  proceed,  we  aust  perform  some  transformations  on  the  IBVP 

(3.1) ,  (3.3),  (3.4).  Through  a  process  described  at  the  beginning  of 
Section  S.2,  the  effects  of  the  initial  data  can  be  incorporated  into  the 
forcing  term  F.  Ve  thus  assuae  f  ■  0.  Ve  also  assume  that  C  «  0  and 
that  the  coefficients  are  constant.  (See  the  comments  after  (4.8))  Ve  now 
apply  a  Fourier  transfora  in  the  tangent  variable  y,  with  dual  variable 

iw,  and  a  Laplace  transfora  in  t,  with  dual  variable  s  =  n  +  i£  (q  >  0). 
The  latter  transform  is  defined  by 

(4.1)  Lw(s)  *  /  e-,tw(t)  dt. 

o 

Under  the  above  transformations  the  systea  (3.1)  becomes 

su(x,a,s)  «  An  +  ?  i<« .BjU  +  F(x.w.s) 

x  j*l  J  J 

Here  we  used  the  fact  f  **  0.  The  IBVP  can  then  be  written  as 
(s)  nx  “  M(a, s)u  -  A_1&  ,  x  >  0 

(4.2) 

(b)  a  I(0.», s)  »  Su*1  +  g(», s). 


where 

(4.3)  M(tt.s)  -  A“1(sl  -  ^  iw-B .) 

j  J  J 

It  aay  not  be  apparent  why  one  would  use  transforms  in  (y,t)  rather 
than,  for  example,  in  (x,y).  This  question  will  be  discussed  at  the 
beginning  of  Section  4.3. 
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For  later  reference,  we  show  that  the  nse  of  Laplace  and  Fonrier 
transforms  enables  one  to  express  solutions  of  (3.1)  as  superpositions  of 
elementary  wave-like  solutions.  The  Laplace  transform  (4.1)  of  a  function 
w  is  the  same  as  the  Fourier  transform  of  the  function 


The  factor  e-,,t  is  inserted  to  ensure  integrability.  The 
application  of  an  inverse  Fourier  transform  shows  that  w(t)  is  a 
superposition  of  modes  e^eH*  *  est ,  where  q  is  fixed. 

The  Fourier  transform  in  y  introduces  an  additional  oscillatory  factor 
exp(iw*y).  Furthermore,  the  solutions  of  the  inhomogeneous  system 
(4.2) (a)  are  based  upon  factors  exp(<x)  for  various  eigenvalues  K  of 
M(u,s).  (See  Section  3.2)  The  solutions  of  the  partial  differential 
equation  are  thus  built  up  from  modes 

(4.4)  eKX  +  iw’y  +  #t  . 

These  modes  will  be  discussed  extensively  later. 

The  normal  modes  (4.4)  were  also  used  in  the  discussion  of  the 
necessary  condition  in  Section  3.  However,  we  deliberately  avoided 
superimposing  such  modes  in  order  to  avoid  unnecessary  complications. 

We  now  describe  the  main  topics  to  be  discussed  in  the  present  section. 

Definition  .  The  characteristic  variety  of  the  system  (3.1)  is  the  set 
of  all  a  e  IR,  <■>  e  IK",  ^  el  such  that 

(4.5)  det[Cl  -  (oA  +  l  (d.B.)]  -  0 

j  j  j 

The  matrix  function  oA  +  1  UjBj  is  commonly  called  the 
"principal  symbol"  of  (3.1).  Some  remarks  in  Section  3.1  imply  that  if 
the  system  is  hyperbolic,  then  the  symbol  has  real  eigenvalues  and  a 


complete  set  of  eigenvectors.  In  general,  A  end  the  Bj  may  depend 
on  (x,y,t),  and  thus  the  characteristic  variety  may  also  depend  on  (x,y,t). 

A  typical  relation  (4.5)  is  graphed  in  Figure  4.1  in  Section  4.2. 

There  are  a  couple  of  reasons  why  the  relation  (4.5)  is  relevant  to  the 
study  of  well-posedness.  First,  a  topic  of  major  interest  in  the  theory  is 
the  structure  of  the  matrix  M(u, s)  in  (4.3).  For  reasons  to  be  discussed 
later,  a  particularly  important  topic  is  the  behavior  of  M(w, s)  as 
Re  s  -*■  0.  This  is  already  suggested  by  the  manner  in  which  the  necessary 
condition  (3.18)  is  strengthened  by  the  D.K.C.  (3.20).  For  the  moment  let 
us  consider  the  limiting  case  M(», i£) .  If  is  it  a  purely  imaginary 
eigenvalue  and  z  is  a  corresponding  eigenvector,  then 

(4.6)  (a)  ioz  *  A”*(i£  -  j[  iWjBj)z, 
or 

(4.6)  (b)  £z  *  (oA  +  £  UjBj )  z. 

The  point  (a,u,£)  must  satisfy  the  relation  (4.5).  The  characteristic 
variety  will  thus  enable  one  to  visualize  the  behavior  of  ll(«»,s)  as 
Re  s  -*-0.  The  relations  (4.6)  may  be  compared  with  the  conclusions  of 
Proposition  3.2. 

A  second  reason  for  studying  the  characteristic  variety  is  that  it  is 
related  to  the  phyaical  phenomenon  of  wave  propagation.  Suppose  that  (3.1) 
has  constant  coefficients  and  that  C  =»  0  and  F  *  0,  so  that  (3.1)  ia 


(4.7)  ut  *  Au_  +  £  B, — 

*  X  jd7j 

If  a  plane  wave  solution 

(4.8)  z  ei(cx  +  «*y  +  £t) 

is  inserted  into  (4.7),  where  z  is  a  vector,  the  result  is  (4.6) (b).  The 
characteristic  variety  thus  describes  the  set  of  all  wave  numbers  and 
frequencies  of  plane  wave  solutions  of  (4.7). 
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This  analysis  is  approximately  valid  in  more  general  eases  if  the 
frequencies  are  sufficiently  large.  When  C  ^  0,  one  obtains  the  equation 

£x  *  («A  +  I  WjBj)z  ♦  Cz 

instead  of  (4.6) (b),  but  the  effects  of  C  are  small  for  large  o.m.C. 

Also,  in  the  osse  of  variable  coefficients  the  coefficients  appear  nearly 
constant  to  a  wave  whose  frequencies  are  sufficiently  high. 

One  may  object  that  the  mode  (4.8)  does  not  satisfy  the  homogeneous 
initial  condition  used  to  obtain  (4.2) (a).  This  issue  is  settled  by  the 
discussion  appearing  in  Section  5.2. 

In  the  theory  of  well-posedness  one  actually  encounters  the  modes  (4.4) 
rather  than  the  modes  (4.8).  However,  the  relations  between  the  two  cases 
are  of  major  interest  and  will  be  discussed  extensively  later.  (See 
Proposition  4.2  and  Section  5.3) 

We  now  outline  the  contents  of  the  remainder  of  this  section.  In 
Section  4.1  we  give  some  examples  of  first-order  hyperbolic  systems.  In 
Section  4.2  we  discuss  the  structure  of  the  characteristic  variety  for 
systems  like  those  mentioned  in  Section  4.1.  In  Section  4.3  we  nse  the 
characteristic  variety  to  describe  some  properties  of  M(tt.i^),  and  we 
point  out  that  the  examples  model  some  crucial  aspects  of  the  behavior  of 
more  general  systems.  A  major  point  of  interest  in  Sections  4.2  and  4.3  is 
the  problem  of  identifying  the  incoming  and  outgoing  portions  of  the 
solution.  In  Section  4.4  we  discuss  the  manner  in  which  these  portions  are 
labelled  algebraically  by  certain  processes  appearing  in  the  theory  of 
well-posedness.  Some  of  the  content  of  Sections  4.2  and  4.3  has  already 
appeared  in  [7]. 

4.1  Pggwp^pi  si  first-order  hyperbolic  systems 

In  this  sub-section  we  list  the  shallow  water  equations,  the 
two-dimensional  isentropic  Euler  equations  of  gas  dynamics,  and  Maxwell's 
equations. 

The  shallow  water  equations  can  be  written  in  the  form 


i 
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(4.9) 


nt  +  oux  +  vuy  +  4X  =  0 

vt  +  nvx  +  vvy  +  4y  -  0 

♦j  +  (40)X  +  (#V)y  “  0 

(For  a  derivation  tee,  e.g.,  Conrant  and  Friedricbs  [2],  Whitham  [31].) 

Thia  systen  describes  tbe  notion  of  an  incoapressible  fluid  when  tbe  fluid 
is  bonnded  below  by  a  solid  boundary  and  its  depth  is  snail  relative  to  a 
typical  horizontal  length  scale.  In  (4.9)  u(x,y,t)  and  v(x,y,t)  are  the 
horizontal  coaponents  of  velocity  and  4  is  tbe  geopotential;  i.e., 

4  *>  gh(x,y,t),  where  g  is  the  acceleration  due  to  gravity  and  h  is  the 
height  of  the  free  snrface.  Coriolis  effects  can  be  included  in  (4.9)  by 
adding  the  terns  -fv  and  fu  to  the  left  sides  of  the  first  and  second 
equations,  respectively. 

The  systen  (4.9)  can  be  written  as 

/»  0  l\  /*  »  °\ 

(4.10)  ,,  ♦  0  .  0  U,  *  0  .  1  ,,-0. 

\  ♦  0  «  /  \  0  4  v/ 

where  q  “  (u,v.4)T.  If  we  pre-nultiply  (4.10)  by  the  diagonal  natri., 
diag  (4.4,1}.  the  result  has  the  fora 

Vt  +  V*  *  Vy  *  °» 

where  the  Aj  are  synnetric  and  A0  is  positive  definite.  (4.10) 
is  therefore  a  " syanetr izable' *  hyperbolic  systea  (e.g.,  Friedrichs  [4], 
John  [9]). 

The  systen  is  nonlinear  and  thus  not  covered  by  the  well-posedness 
theory  which  is  the  subject  of  the  present  paper.  However,  the  theory  is 
applicable  to  the  linearised  version  of  the  systen.  To  obtain  this  version, 
suppose  that  q  *  (u,v,4)T  and  q  ♦  q'  -  (a  +  u',  v  ♦  v'.  4  +  4')T 
are  solutions  of  (4.10),  and  suppose  that  the  perturbations  u',v’,4* 
are  snail.  If  we  substitute  q  q*  into  (4.10)  and  neglect  terns 
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which  arc  quadratic  in  the  perturbation,  the  result  is  the  following  linear 
system  for  q' : 


(4.11)  q't  + 


0 


v 

4 


-  o 


The  matrix  D  involves  various  derivatives  of  n.v.p;  and  if  the  system 
(4.9)  includes  Coriolis  effects  then  the  parameter  f  would  also  appear  in 
D. 

Ve  now  transform  (4.11)  to  a  symmetric  system.  Ve  first  note  that 
(4.11)  is  symmetrizable  in  the  same  sense  that  the  nonlinear  system  (4.10) 
is  symmetrizable.  This  would  yield  the  form  (3.1)  discussed  earlier,  with 
A  and  the  Bj  symmetric,  except  that  the  time  derivative  would  have  a 
coefficient  matrix  which  is  not  the  identity.  In  order  to  avoid  this 
coefficient,  we  introduce  a  change  of  variable.  Let  Q  =  diag  {/$,  ^4,  1) , 
and  note  that 

(u  0  1  \  /  u  0  \ 

0  u  0  Q_1  -  I  0  u  0 

4  0  u  /  \  /4  0  u  / 


A  similar  conclusion  can  be  reached  about  the  coefficient  of  q’^.  The 
system  (4.11)  can  therefore  be  transformed  to 


(4.12)  wt  + 


(u  0  /4  \  /  v  0  0  \ 

0  u  0  J  wx  +  Jo  v^Jwy  +  Ew 

^  0  u  /  \  0  /4  v  / 


0. 


where  w  =  Qq*  *  ( u '  /(T,  v'^,  4').  The  coefficient  E  is  different 
from  the  coefficient  D  in  (4.11)  because  of  some  derivatives  of  Q  which 
are  encountered  during  the  change  of  variable.  The  quantity  ^  *  /gh"  is 


the  speed  of  gravity  waves  in  shallow  water  and  will  be  denoted  below  by  c. 

The  two-diaiensional  isentropie  Enter  equations  of  gas  dynasties  have  a 
structure  which  is  similar  to  that  of  the  shallow  water  equations.  The 
Euler  equations  are 

1  dp 

ut  +  uu_  +  vn  *  - 

1  y  p  dx 

1  dp 

(4.13)  vt  +  nv_  +  wv  - - 

y  p  3y 

Pt  +  (pu)x  +  (pu)y  -  0 
p  "  p(p) 

Here  u(x.y.t)  and  v(x.y.t)  are  the  x-  and  y-components  of  velocity', 
respectively;  p  is  the  density  of  the  fluid;  and  p  is  the 
pressure.  (4.13)  can  also  be  written  as 


(4.14) 


TP 

where  q  *=  (n.v.p)  .  The  systesi  (4.14)  can  be  symmetrized  by 

1  dp  dp 

pre-mult  ip lying  by  the  diagonal  matrix  diag  (p.p.  ~  if>) .  If  dp  >  0 

this  siatriz  is  positive  definite,  so  that  (4.14)  is  symmetrizable 
.  ...  ^2  .  _  .... 


hyperbolic.  From  now  on  we  will  assume 


>  0;  in  many  situations 


p  »  KpY.  where  y  and  K  are  constants  and  y  >  1. 

The  system  (4.14)  can  be  linearized  and  then  symmetrized  to  produce  a 
form  like  that  of  (4.12).  This  process  yields 


(4.15) 


0  0 


wy  +  Ew  =  0 


C 

Here  w  ■  (u'.v*.  (-)p’)  ,  where  u*  and  v'  are  the  perturbations  in 


the  velocity  components,  p*  is  the  perturbation  in  the  density,  and 
/3p 

c  m  *  "Sp  •  The  latter  quantity  is  the  loeal  speed  of  sound. 

The  above  disonssion  implies  that  the  linearised  shallow  water 
equations  and  the  linearized  Euler  equations  can  each  be  written  in  the  fora 
(3.1)  discussed  earlier, 

(4.16)  wt  «  Awx  +  BWy  +  Cw, 

where  in  this  case 


and  C  ■  -E.  An  analogous  form  can  be  obtained  for  the  Euler  equations  in 
three  dimensions. 

For  later  reference  we  discuss  briefly  the  relation  (4.5)  for  the 
system  (4.16).  The  principal  symbol  of  (4.16)  is  oA  +  wB,  and  its 
eigenvalues  are 


(4.17) 


■  -(u,v)*(o.w)  *  -uo  -  v» 

^2*  ^3  *  -(u.v)*(o,»)  +  c(c2  +  ((2)1/2 


A  typical  relation  (4.17)  is  graphed  in  Figure  4.1.  The  eigenvectors 
corresponding  to  the  eigenvalues  are 


(4.18) 


*  (-«,o,0)T 


z2,z3  “  (±»»  ±“.  <«2  +  w2)1^2)* 


Ve  now  discuss  Maxwell's  equations  of  electromagnetism.  These  are 


(4.19) 


(a) 

(b) 

(c) 


div  E  «  — 


div  B 


dE 


—  ■  c^(curl  B  -  p.J) 
3t  0 


(d) 


9B 

—  ■  -  curl  E, 

at 


where  E  is  the  electric  field  intensity,  B  is  the  magnetic  induction, 
p  is  the  total  charge  density,  J  is  the  current  density,  c  is  the  speed 
of  light,  and  s0  end  p0  are  constants.  E,B  and  J  are  vectors 
having  three  components. 

In  certain  situations  the  system  (4.19)  can  be  simplified  by  deleting 
the  first  two  equations.  The  reason  is  the  following  (e.g.,  Lorrain  and 


Corson  (121).  If  we  take  the  divergence  of  (d),  the  result  is 

a 

(div  B)  -  -  div  curl  E  *  0.  If  at  each  point  in  space  div  B  ■*  0 


at  some  time,  then  (b)  follows  from  (d)  and  may  be  omitted.  Similarly, 

equation  (c)  implies  (div  E)  *  ~«^M0  div  J  *  c^p0  dp/dt  *  dp/3t, 

—1/2  P 

since  c  «  (e  p  )  •  Tb®*  div  E  and  —  differ  by  a  constant  at 

o 

each  point  in  space.  If  for  each  point  they  are  simultaneously  zero  at  some 


i 


time,  then  equation  (a)  is  a  consequence  of  (c)  and  may  also  be  eliminated  from 
the  system. 

Under  the  above  conditions  Maxwell's  equations  reduce  to  the  form 


r 


It  follows  that  (4.20)  it  symmetric  hyperbolic  if  c  =  1.  If  c  4  1  the 
symmetric  fora  cen  be  obtained  by  replacing  B  with  c_1E  and  t  with 
ct . 

The  principal  syabol  of  (4.20)  will  be  of  interest  in  later 
discussions.  A  calculation  shows  that  its  eigenvalues  all  have  multiplicity 
two  and  are  given  by 

^1'  ^2  “  0 

(4.21)  K3.  C4  -  c(o2  +  Ul2)1/2 

£5.  £6  «  -c(o2  +  lewl2) 1/2 

1  i2  2  2 

Here  ,re  dual  to  re*Pect ively,  and  lul  ”  uij  +  w2 . 

The  coaponents  of  the  corresponding  eigenvectors  are  rational  functions 
of  the  quantities  o.Wj.Wj  and  (o2  +  l«l2)2^2. 


4.2  Structure  of  the  characteristic  variety;  troun  velocity. 

In  this  section  we  discuss  the  structure  of  the  characteristic  variety 
and  how  this  structure  is  related  to  the  problem  of  identifying  the  incoming 
and  outgoing  portions  of  the  solution.  As  suggested  in  Section  1,  the 
latter  question  can  be  of  fundamental  importance  in  the  study  of 
well-posedness  of  initial-boundary  value  problems. 

A  typical  relation  (4.5)  is  graphed  in  Figure  4.1.  This  is  the 
characteristic  variety  for  the  linearized  shallow  water  equations  and  the 
linearized  two-dimensional  Euler  equations.  (See  (4.17))  The  plane 
corresponds  to  the  eigenvalue  £2,  and  the  cones  correspond  to  £j 
and  £3.  The  cones  are  right  circular  cones  if  and  only  if  the  plane  is 
horizontal.  We  let  ft  denote  the  double  cone  and  T  denote  the  set  of  all 
(«.£)  which  correspond  to  points  on  ft . 


(a)  Characteristic  variety  (b)  Projection  onto  («,4 )  space 

Figure  4.1 

The  characteriatic  variety  for  Maxwell's  equations  (4.20)  has  a  similar 
structure.  In  this  case  each  detensination  of  4  is  a  function  of  three 
variables  (a,*^,*^),  and  each  snch  4  has  sinltiplicity  two. 

The  ssiallest  eigenvalue  is  sero  and  would  thus  correspond  to  a  horizontal 
plane  in  Figure  4.1(a). 

Aa  noted  in  the  introduction  to  Section  4,  the  relation  (4.5)  describes 
the  set  of  all  (e,n,4)  in  wave-like  solutions 

(4.22)  #i(ox  +  t»*y  ♦  4t) 

of  the  systen  (4.7),  which  is  a  special  case  of  (3.1).  In  the  current 
discussion  we  are  particularly  interested  in  the  directions  in  which  these 
waves  propagate.  In  order  to  deal  with  this  question,  we  first  make  some 
renarks  about  the  concept  of  group  velocity  as  it  applies  to  the  present 
situation.  We  then  show  how  incoming  and  outgoing  portions  of  the  solution 
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can  be  identified  directly  froa  graphs  like  Fignre  4.1(a). 

This  discussion  will  be  a  description  of  solutions  of  (3.1)  when  C  *  0 
and  the  coefficients  are  constant  ;  and  for  sufficiently  high  frequencies 
it  will  be  an  approxinate  description  of  the  general  ease,  where  C  aay  be 
nonzero  and  the  coefficients  aay  depend  on  (x,y,t).  (See  the  introduction 
to  Section  4.)  At  the  end  of  this  sub-section  we  relate  this  discussion  to 
the  theory  of  propagation  of  singularities,  which  provides  a  rigorous 
treataent  of  the  aore  general  case. 

For  each  fixed  t,  the  individual  aode  (4.22)  is  constant  along  lines 
in  the  (x.y)  plane  for  which  (x,y)*(o,a)  is  constant.  It  then 
follows  that  the  wave  propagates  in  the  direction  (o.a)  with  phase  speed 
~5(o2  +  M2r1/2.  However,  when  one  superiaposes  various  waves  of  the 
fora  (4.22),  the  net  effect  is  to  produce  a  solution  in  which  energy  is 
propagated  with  a  velocity  which  aay  be  different  froa  the  phase  velocities 
of  the  individual  waves.  This  velocity  is  the  group  velocity,  and  in  the 
notation  of  (4.22)  it  is  given  by 
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(4.23)  grad(-£)  -  ( - - ) 

do  3« 

Group  velocity  and  phase  velocity  differ,  for  exaaple,  in  wave  groups 
which  propagate  in  one  diaension  and  which  are  "dispersive”,  i.e.,  the 
phase  velocity  varies  with  wavelength,  (e.g.,  Vhithaa  [31])  In 
mult-i-diaensional  problems  the  two  velocities  aay  also  differ  for  reasons 
which  are  not  related  to  variations  in  wavelength.  This  effect  can  be  seen 
aost  easily  in  the  case  of  the  siaple  equation 

(4.24)  ut  +  ux  +  Uy  =  0 

The  solutions  to  this  equation  are  constant  along  characteristic  corves  for 
dx  dy 

which  <fF  “  TF  *  *nd  tbey  thus  represent  translations  in  the  (1,1)  direction. 
However,  when  one  substitutes  a  wave  fora  (4,22)  into  (4.24),  the  result 
is  C  *  —a—u,  and  by  choosing  (o,«)  appropriately  one  can  obtain  a  wave-like 
solution  having  a  phase  velocity  which  points  in  any  desired  direction. 

This  direction  cannot  be  aeaningful  because  of  the  translatory  nature  of 

(4.24) .  In  fact,  suppose  that  a  family  of  such  waves  is  superimposed  by  the 


integral 


/'Hex  +  «y  -  <a+«)t)  f(a  u)  dadUf 

■A 

where  f  is  an  amplitude  function.  This  integral  can  be  written 
Jei[  o(x-t)  +  «(y-t) ]  f(a  u)  d(f  du 

-  f  (x  -  t,  y  -  t) , 

where  f  is  the  inverse  Fourier  transform  of  f.  The  individual  waves  can 
■ove  with  phase  velocities  having  any  direction  whatever,  but  the  waves 

superpose  to  produce  a  simple  translatory  effect.  In  this  case  the  group 

3  fj  9E, 

velocity  is  (-  jj,  -  *  (1,1),  which  is  the  direction  of  translation. 

Some  interesting  applications  of  the  concept  of  group  velocity  are 
found  in  the  study  of  finite  difference  approximations.  See,  e.g., 
Trefethen  [27],  [28]  and  Vichnevetsky  and  Bowles  [30]. 

Ve  now  relate  the  directions  of  the  group  velocity  to  the  structure  of 
the  characteristic  variety  graphed  in  Figure  4.1(a).  In  this  particular 
example  there  are  three  determinations  of  £  as  a  function  of  (o.u). 

The  determination  corresponding  to  the  plane  is  simplest  to  analyze;  its 
gradient  is  constant,  so  the  group  velocity  is  a  constant  function  of 
(o,«).  The  plane  thus  represents  a  translatory  motion  which  is  similar 
to  that  found  in  solutions  of  the  equation  ut  +  ux  +  Uy  =  0. 

In  the  example  illustrated  here  the  group  velocity  points  into  the  spatial 
domain  x  >  0. 

The  cones  are  more  complicated  than  the  plane.  In  Figure  4.2  we 

show  cross-sections  of  these  cones  for  fixed  £  >  0  and  fixed  £  <  0.  The 

directions  of  group  velocity  are  indicated  by  the  solid  arrows;  the  group 

JC  3C 

velocity  is  grad  (-5)  *  (-  ^j) ,  so  the  direction  of  this  velocity 

is  opposite  the  direction  of  most  rapid  increase  in  £  «=  C(o,w).  These 

directions  are  orthogonal  to  the  level  sets  shown  in  the  Figure.  In  this 

picture  the  group  velocities  are  different  from  the  phase  velocities,  since 
the  latter  always  point  along  lines  which  pass  through  the  origin. 
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In  parts  (a)  and  (b)  of  Figure  4.2  the  points  A,B  denote  the  points 
(«»  <o,  £ )  corresponding  to  group  velocity  which  is  tangent  to  the 
boundary  z  ■  0.  The  relevance  of  these  points  will  be  discussed  below  and 
in  the  next  sub-section.  Ve  now  explain  their  locations.  The  vertical 
coordinates  of  points  on  ft  are  given  by  susis  of  vertical  coordinates  of 
the  plane  and  those  of  a  right  circular  cone.  (See  (4.17))  At  each  point 
A  or  B  the  gradient  of  the  right  circular  cone  is  denoted  by  the  dotted 
arrow  C,  which  points  along  a  line  passing  through  the  origin.  The 
gradient  of  the  plane  is  indicated  by  the  dotted  arrow  P.  Their  suai  lies 
along  the  sane  line  as  the  group  velocity. 

The  cones  ft  correspond  to  the  propagation  of  sound  waves  and  gravity 
waves,  respectively,  for  the  linearized  Enler  equations  and  the  linearized 
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•hallow  water  equations.  In  the  analogous  picture  for  Maxwell's  equations  (2 
would  correspond  to  the  propagation  of  electromagnetic  waves. 

We  now  mention  some  algebraic  aspects  of  the  discussion  given  above. 

We  have  denoted  by  T  the  set  of  all  («,£)  corresponding  to  points  on  & 
i.e..  r  is  the  projection  of  &  onto  the  (», £)  space.  For  the  example 
illustrated  here  the  cones  ft  determine  two  values  of  o  for  each  (»,  £) 
in  the  interior  of  T.  The  above  discussion  implies  that  each 
determination  of  a  m  o(w.C)  on  each  component  of  F  can  be 
associated  with  motion  into  or  out  of  the  spatial  domain  x  >  S. 

Furthermore,  some  remarks  made  in  the  introduction  to  Section  4  imply  that 
for  each  such  «.  the  number  io  is  an  eigenvalue  of  M(t»,  i£).  (See 
(4.6) (a))  This  suggests  that  one  can  identify  incoming  and  outgoing 
portions  of  the  solution  by  diagonalizing  M(u,iO.  This  subject  will  be 
addressed  in  the  next  sub-section.  As  part  of  this  discussion  we  point  out 
that  M(w.  iC)  is  defective  when  («s.  £)  corresponds  to  tangential 
incidence.  In  Section  4.4  we  describe  a  process  by  which  the  incoming  and 
outgoing  portions  of  the  solution  are  labelled  algebraically  in  the  theory 
of  well-posedness. 

We  now  relate  the  earlier  discussion  of  group  velocity  to  the  theory  of 
propagation  of  singularities.  The  latter  is  mainly  a  discussion  of  the 
propagation  of  high-frequency  portions  of  the  solution,  since  the  smoothness 
of  a-  function  is  governed  by  the  rate  of  decay  of  its  Fonrier  transform. 

Denote  (x.y.t)  by  z  and  (o.m.  £)  by  £  ,  and  let  p(z.C)  denote  the 
determinant  in  (4.5).  The  singularities  in  the  solution  propagate  along 
bicharacteristic  curves,  which  are  curves  (z(t).C(t))  satisfying  the 
Hamilton-Jacobi  equations 

(a)  z  *  gradj.  p(z,C) 

(4.25)  (b)  1  -  -gradz  p(x.C) 

and  the  constraint  p(z.O  “  0.  (e.g.,  Nirenberg  [17],  Taylor  [26])  The 

curve  z(x)  denotes  the  path  of  propagation  in  the  (x.y.t)  space,  and  the 
curve  5(t)  records  any  frequency  shifts  which  take  place  during  the 
propagation.  For  a  discussion  of  ray  theory  in  a  more  applied  setting,  see, 
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e.g.,  Whit ham  [31]. 

Equation  (4.25) (a)  ttys  that  the  direction  z  ia  parallel  to 
grad^  p(z,C),  i.e.,  ia  orthogonal  to  the  level  aet  (£:p(z,£)  11  0).  This 
set  is  the  characteristic  variety  defined  by  (4.5);  an  example  is  graphed 
in  Figure  4.1.  It  then  follows  that  the  projection  of  z  onto  the  spatial 
variables  (x,y)  is  perpendicular  to  level  sets  of  the  type  graphed  in  Figure 
4.2.  That  is.  the  direction  of  propagation  liea  along  the  Bane  line  as  the 
group  velocity.  One  can  check  easily  that  the  direction  of  propagation  is 
the  direction  of  the  group  velocity,  rather  than  the  exact  opposite;  for 
exanple,  on  the  upper  cone  (  £  >0)  in  Figure  4.1,  an  increase  in  t 
corresponds  to  a  movement  of  (x.y)  toward  the  center  of  the  cone. 

4.3  Properties  of  M(w,iC) 

In  this  sub-section  we  use  the  structure  of  the  characteristic  variety 
to  describe  some  properties  of  the  matrix  M(w,iO  ■=  A“*(i£I  -  £  iWjBj). 

(See  (4.3))  Ve  will  deal  mainly  with  the  diagonalizability  of  this  matrix  and 
the  occurrence  of  purely  imaginary  eigenvalues.  A  noted  in  the  previous 
sub-section,  the  behavior  of  *!(«, i£)  ia  closely  related  to  the  problem  of 
identifying  the  incoming  and  outgoing  portions  of  the  solution. 

The  matrix  M(w, s)  arises  from  the  use  of  Fourier  and  Laplace 
transforms  with  respect  to  y  and  t,  respectively.  Before  we  proceed  with 
the  discussion  of  M(u, i£),  we  first  mention  why  one  would  transform  in 
(y.t)  rather  than  in  some  other  set  of  variables. 

First,  there  can  be  a  problem  with  transforming  in  the  normal  variable 
x.  Such  a  transformation  would  require  information  about  the  solution  away 
from  the  boundary  x  **  0,  and  this  does  not  seem  appropriate  in  a  discussion 
of  boundary  conditions.  This  difficulty  can  be  avoided  by  transforming  in 
(y.t) . 

A  second  reason  is  related  to  the  problem  of  identifying  incoming  and 
outgoing  modes.  Ve  have  already  noted  that  the  eigenvalues  io  -  io(«*,£) 
of  M(w, i£)  are  very  useful  in  this  regard.  On  the  other  hand,  if 
we  were  to  transform  (3.1)  with  respect  to  (x,y),  say,  then  the  result  would 
be 
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(if  C  =  0,  F  «*  0) .  The  properties  of  solutions  of  this  system  are  governed 
by  the  eigenvalnes  i£  of  ierA  +  £iu>.Bj.  Unfortunately,  the 
discussion  in  Section  4.2  implies  that  various  incoming  and  outgoing  modes 
are  mixed  together  in  certain  determinations  of  i£  as  a  function  of 
(o.w).  These  determinations  are  associated  vith  the  cones  ft.  It  thus 
does  not  seem  possible  to  accomplish  the  desired  separation  by  studying  the 
transformed  equation  (4.26). 

We  now  discuss  the  eigenvalues  of  M(te,iO.  We  recall  that  if  io 
is  a  purely  imaginary  eigenvalue  and  z  is  a  corresponding  eigenvector, 
then 

(4.27)  iorx  =  A_1(ia  -  }jiWjBj)z, 
or 

(4.28)  i£z  **  (ieA  +  £is>jBj)z 

According  to  remarks  made  in  the  introduction  to  Section  4,  the  symbol 

(4.29)  ioA  +  £ iw.Bi 
j  j  j 

has  purely  imaginary  eigenvalues  and  a  complete  set  of  real  eigenvectors. 
This  fact  will  be  exploited  in  Proposition  4.1. 

In  general,  the  number  of  purely  imaginary  eigenvalues  of  M(w.i^) 
may  vary  with  the  position  of  («,£).  In  the  case  of  the  shallow  water  or 
Euler  equations  there  are  two  eigenvalues  io  associated  with  the  double 
cone  Q,  when  (w, Q  is  in  the  interior  of  r .  (See  Figures  4.1  and  4.2) 

As  (w, O  approaches  the  edge  of  T,  these  eigenvalues  coalesce,  and  as 
(«,£ )  leaves  T  these  eigenvalues  leove  the  imaginary  axis.  They  cannot 
be  purely  imaginary,  since  for  any  imaginary  eigenvalue  io  the  point 
(o, )  must  lie  on  one  of  the  surfaces  in  Figure  4.1. 

In  fact,  for  («>,£)  outside  T  these  eigenvalues  must  have  the  form 
±  p  +  io,  where  p  is  real.  The  eigenvalues  of  the  real  matrix 


i 0  -  A-1(Cl  -  IwjBj) 

are  either  real  or  cone  in  complex  conjngate  pairs  o  +  ip.  Vhen  we 
multiply  by  1  the  result  is  +p  +  la. 

For  the  Euler  and  shallow  water  equations  there  is  also  a  value  io 
associated  with  the  plane  illustrated  in  Figure  4.1.  In  many  cases  there  is 
no  difficulty  with  expressing  a  as  a  function  of  («,€) .  The  plane 
satisfies  the  equation 

£  *  -  on  -  u>v 

(see  (4.17)),  so  ia  can  be  written  in  terms  of  (u,0  whenever  u  ^  0. 

This  condition  is  equivalent  to  requiring  det  AM.  where  A  is  the 
coefficient  of  wx  in  (4.16).  In  Section  3.1  we  assumed  that  this 
condition  was  satisfied,  and  we  used  it  to  write  M(w,s)  *  A_*(sl  -  I  iWjBj) 
(e.g.,  (3.9)). 

The  assumption  det  A  /  0  has  a  physical  interpretation.  The  vector 
field  (u,v)  associated  with  the  system  (4.16)  is  the  velocity  of  the 
flow  about  which  the  system  has  been  linearised.  If  A  is  to  be 
nonsingular,  then  this  flow  cannot  be  zero  and  cannot  be  tangent  to  the 
boundary  x  =  0.  In  Figure  4.1(a)  the  plane  cannot  be  horizontal  and  cannot 
have  a  gradient  which  is  parallel  to  the  w-axis.  The  assumption  det  A  t  0 
is  thus  not  always  valid  for  the  shallow  water  and  Euler  equations. 

In  addition,  it  is  never  valid  for  Harwell's  equations.  This  can  be 
seen  from  the  representation  (4.20)  for  the  system  or  from  the  formula 
(4.21). 

Ve  now  make  some  remarks  about  the  diagonalizabil ity  of  VHu.iQ. 
Suppose  that  the  characteristic  variety  has  a  cone  ft  ,  and  suppose  that  for 
(«>,  Q  in  T  there  is  no  difficulty  with  solving  for  values  of  io 
associated  with  surfaces  other  than  ft.  There  will  be  no  need  to  make  any 
assumptions  regarding  the  multiplicities  of  any  of  the  eigenvalues.  The 
following  Proposition  may  be  contrasted  with  Proposition  3.2,  which 
describes  M(«i,s)  for  Re  s  >  0. 
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Pronos  it  Ion  4.1.  If  (»,£)  Is  in  the  interior  of  T  ,  then  M(t»,iO 
=  A“*(i£I  -  ^ietjBj)  has  purely  inaginary  eigenvalues  and  a 
complete  set  of  real  eigenvectors.  This  is  not  the  case  if  (»,£> )  is 
outside  F  .  The  eigenvectors  can  be  determined  from  those  of  the  symbol 

(4.29) .  ioA  + 

Proof.  Equations  (4.27)  and  (4.28)  show  that  the  eigenvectors  of 

(4.29)  are  also  eigenvectors  of  M(u,iC).  We  know  that  (4.29)  has  a 
complete  set  of  real  eigenvectors  corresponding  to  fixed  (o.«o)  and 
various  eigenvalues  i£.  We  want  to  shov  the  same  thing  for  M(w, iO«  f or 
fixed  (u,0  in  T  and  various  eigenvalues  ic. 

Suppose  that  <«,£)  is  in  T  .  and  let  Oj,...  denote 
the  eigenvalues  of  M(w,  i£) .  For  each  Oj  choose  a  basis  Ej 
for  the  eigenspace  of  iojA  +  JiwjBj  corresponding  to  the 
eigenvalue  i£.  We  are  allowing  for  the  possiblity  that  (4.29)  might  have 
multiple  eigenvalues.  The  elements  of  Ej  are  also  eigenvectors  of 
M(w,H)  corresponding  to  the  eigenvalue  ier .  We  claim  that  the  union 
of  the  Ej  is  a  complete  set  of  vectors.  There  are  clearly  enough  of 
these  vectors.  The  fact  that  they  are  linearly  independent  follows  from  an 
argument  which  is  essentially  the  one  which  shows  that  eigenvectors 
corresponding  to  distinct  eigenvalues  axe  linearly  independent.  This 
completes  the  proof. 

We  make  some  comments  about  the  behavior  of  M(w, i  £)  when  («,£) 
lies  on  the  edge  of  T.  According  to  Figures  4.1  and  4.2,  this  case 
corresponds  to  group  velocity  which  is  tangent  to  the  boundary,  and  it  also 
corresponds  to  the  coalescence  of  different  values  of  o( »,£).  In 
addition,  as  (*>,  0  approaches  the  edge  of  T  ,  various  eigenvectors  of 
M(<i>,  i£)  associated  with  the  two  determinations  of  io(w.C)  come 
together,  so  that  N(u,iQ  fails  to  have  a  complete  set  of  eigenvectors. 

In  the  case  of  the  Euler  and  shallow  water  equations,  M(«, i£)  acquires  a 
2x2  nondiagonal izable  block.  In  systems  for  which  ft  corresponds  to 
multiple  eigenvalues,  the  nondiagonalizable  block  can  be  larger.  This 
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occurs,  for  example,  with  Maxwell's  equation.  (In  this  case  it  is  not 
possible  to  write  M(t>,s),  since  det  A  ■  0,  but  there  is  an  analogous 
■atrix  which  is  used.  See  Section  1  of  Majda  and  Osher  [14].)  The 
defective  behavior  of  M(u,i3  causes  major  difficulties  in  the  theory  of 
well-posedness.  (See  Section  6) 

Proposition  4.1  snd  the  accompanying  remarks  have  obvious  extensions 
to  more  general  cases,  e.g.,  where  the  characteristic  variety  has  several 
different  cones.  In  any  system  the  msjor  point  of  interest  is  the  behavior 
of  M(», iO  when  the  eigenvalues  io  coalesce,  that  is,  at  tangential 
incidence.  The  examples  mentioned  here  contain  this  principal  difficulty. 
For  the  sake  of  definiteness  and  clarity,  we  orient  many  of  the  discussions 
in  the  remainder  of  the  paper  to  these  examples. 

4.4.  An  gjggfrgqjq  llfcglliSl  Si  incoming  gnd  outgoing  Bgfl&j 

In  this  section  we  discuss  a  process  by  which  incoming  and  outgoing 
modes  are  labelled  algebraically  in  the  theory  of  well-posedness.  We 
consider  the  elementary  solutions 

(4.30)  eiffX  +  ltt*y  +  i^t 

of  (4.7)  which  have  been  discussed  extensively  in  previous  sub-sections.  We 
analyze  what  happens  when  the  dual  variables  i£  and  io  are  perturbed 
to  the  complex  values  s  *  rj  +  i£  and  k  =  p  +  io,  respectively,  so 
that  the  resulting  form 

(4.31)  e<x  +  *•**  +  Mt 

is  also  a  solution  of  (4.7).  The  sign  of  Re  k,  when  Re  s  >  0,  will 
indicate  whether  the  mode  (4.30)  corresponds  to  group  velocity  pointing  into 
or  out  of  the  spatisl  domain  x  >  0.  This  labelling  is  used  in  a 
fundamental  manner  in  the  theory  of  well-posedness.  (See  Section  5)  The 
relationships  between  the  modes  (4.30)  and  (4.31)  will  be  discussed  more 
extensively  in  Section  5.3. 


» 
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The  labelling  process  depends  on  an  assumption  that  <  and  s  are 
invertible  analytic  functions  of  each  other  for  q  near  zero.  This 
analytic  dependence  is  found  in  the  examples  which  have  been  discussed 
earlier.  (See  (4.17)  and  (4.21))  However,  it  fails  for  points  (o.»,  £  ) 
on  the  cones  ft  when  («, £  )  is  on  the  edge  of  ^  (i.e.,  tangential 

incidence).  It  also  fails  for  points  corresponding  to  the  plane  in  Figure 
4.1  when  the  plane  is  horizontal.  This  situation  is  encountered  when  the 
boundary  is  characteristic. 

Proposition  4.2.  Suppose  that  k  and  s  have  the  analytic 
dependence  mentioned  above,  and  suppose  that  (4.30)  is  perturbed  so  that 
Re  s  >  0.  If  (4.30)  corresponds  to  group  velocity  pointing  into  the  spatial 
domain  x  >  0,  then  <  is  perturbed  so  that  Re  K  <0.  If  the  group 
velocity  points  out  of  the  domain,  then  Re  <  >0. 

*£ 

Proof .  The  vector  group  velocity  is  (-  3^,  -Vw£)(*ee  (4.23)). 
Incoming  and  outgoing  modes  thus  correspond  to  d£/do  negative  and 
positive,  respectively.  The  former  case  is  illustrated  in  Figure  4.3;  the 
solid  arrows  indicate  that  perturbations  in  £  and  o  have  opposite  signs. 


s-plane 


k -plane 


i i 

Figure  4.3.  The  case  g0  <  0. 


Because  of  the  analytic  dependence,  the  complex  derivative 
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4^  exists  end  is  equal  to  q9.  Perturbations  in  tc  and  s  thus  have 

a  ratio  which  is  independent  of  direction.  In  the  case  of  Figure  4.3  these 

perturbations  are  negatives  of  each  other,  so  that  Re  <  <0  when  Re  s  >  0. 

This  is  illustrated  by  the  dotted  arrows.  Similarly,  if  (4.30)  corresponds  to 

££  di 

an  outgoing  group  velocity,  then  .j  >  0.  In  this  case  >  0,  and 
Re  s  >  0  implies  Re  <  >0.  This  completes  the  proof. 

The  above  proof  resembles  an  argument  given  by  Taylor  [26, p  202]  in  a 
discussion  of  reflection  of  singularities.  A  similar  argument  was  also 
given  by  Trefethen  in  a  study  of  the  relationship  between  group  velocity  and 
the  stability  of  boundary  conditions  for  finite  difference  approximations 
(e.g..  [27],  [29]). 

We  should  note  that  there  exist  modes  (4.31)  which  cannot  be  obtained 
by  perturbing  (4.30).  When  Re  s  4  0  these  approach  the  form  exp(<x  + 
im'y  +  i£t),  where  Re  <  j *  0.  Such  modes  are  found  for  (*>,  Q 
outside  T.  (See  Section  5.3) 
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5.  Interpret** tons  of  the  Uniform  Ireiu  Condition 

In  this  section  we  discnss  some  effects  of  the  ''Unifora  Kreiss 
Condition"  (U.K.C.).  This  is  the  condition  (3.20)  which  wss  introduced  in 
Section  3.3  snd  which  is  sufficient  to  sssure  the  well-posedness  of  the  IBVP 
(3 .1) , (3 .3) . (3 .4) .  In  this  section  wc  show  thst  the  U.K.C.  csn  be  regsrded 
ss  s  solvability  condition;  it  ensbles  one  to  solve  for  certsin 
"incoming"  dependent  vsrisbles  in  teras  of  "outgoing"  vsrisbles  snd 
boundary  data.  The  U.K.C.  is  used  for  this  purpose  in  proofs  of 
well-posedness.  The  above  interpretation  is  an  analogue  to  the  situation  in 
one  space  dimension. 

The  present  discussion  is  based  on  the  use  of  Fourier  transforms  and 
Laplace  transforas.  The  conclusions  reached  here  are  therefore  limited  to 
systems  which  have  constant  coefficients.  Analogous  conclusions  may  be 
expected  for  the  high-frequency  portions  of  solutions  to 

variable-coefficient  problems,  since  the  coefficients  appear  nearly  constant 
to  waves  whose  frequencies  are  sufficiently  high.  This  principle  is 
contained  in  the  theory  of  propagation  of  singularities  (e.g.,  Nirenberg 
[17].  Taylor  [26])  and  in  discussions  of  "slowly  varying"  wavetrains  in 
the  applied  literature  (e.g.,  Whitham  [31]).  However,  the 
variable-coefficient  case  is  technically  more  complicated,  and  it  may  be 
particularly  difficult  near  tangential  incidence.  (See  Section  5.3) 

In  any  case,  the  properties  of  the  high-frequency  portions  of  the 
solution  are  of  major  interest;  estimates  involving  L  norms  and 
Sobolev  norms  (e.g.,  (3.5))  play  a  major  role  in  the  theory,  and  these  norms 
are  governed  by  behavior  at  high  frequencies. 

In  Section  5.1  we  give  an  outline  of  some  effects  of  the  U.K.C.  This 
discussion  is  based  on  the  system  (4.2) (a)  of  ordinary  differential 
equations  which  was  obtained  through  the  use  of  Fourier  and  Laplace 
transforms.  The  structure  of  the  solutions  of  this  system  is  discussed  in 
Section  5.2. 

In  Section  5.3  we  discuss  in  detail  the  nature  of  the  ''incoming  "  and 
"outgoing"  modes  and  how  these  modes  are  affected  by  the  U.K.C.  The  case 
of  tangential  incidence  is  included  in  this  discussion.  The  development  in 


this  sob-section  is  given  in  tens  of  the  shallow  water  equations  and  the 
two-dimensional  Enler  equations. 

In  Section  3.4  we  discuss  soae  weak  fons  of  well-posedneas  which  can 
be  encountered  if  the  U.K.C.  is  not  satisfied  completely. 

5.1.  General  mutl 

In  this  sub-section  we  outline  some  effects  of  the  D.I.C.  The 
discussion  is  based  on  the  transfoned  problem  (4.2), 

(a)  ux(x,w,s)  *  s) u  -  A-1F(x,«,s) 

(5.1) 

(b)  aI(0,w,s)  *  Su11  +  g(w.s), 
which  was  derived  earlier. 

In  order  to  study  the  solutions  of  this  problem  we  transfon  N  to 
block  fon.  Let  Q(«. s)  be  a  matrix  such  that 

/  "l  o 

(5.2)  <TXMQ  -  I 

\  0  lij 

The  transfonation  is  chosen  so  that  the  eigenvalues  ic  of  and 
1*2  satisfy  Re  <  <  0  and  Re  <  >0,  respectively.  The  dimensions 

of  and  are  thus  H  x  H  and  (n  -  4  )  x  (n  -  l  ) .  (cf. 

Proposition  3.2) 

In  the  theory  of  well-posedness  a  great  deal  of  attention  is  paid  to 
finding  block  forms  which  are  smooth  functions  of  (t»,  s).  Ve  will  not 
worry  about  this  now.  In  (5.2)  the  matrices  Q.M^,  ,nd  *2  »*e  not 

determined  uniquely;  this  will  be  discussed  below. 

The  system  (5.1)(a)  can  be  written  as 

Q_1ux  -  Q-1)KHQ_1u  -  Q_1A_1F, 
or 
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(5.3)  vz(x,w,s) 


v  +  G, 


“2 


Q  *u.  For  convenience  we  partition  v  into  vectors  vJ 


.v,)J  and 


v11  =  <v0 


»  •  •  •  '  • 


,vn)  can  be 


where 

-  <▼! . ,  .»»  >  - 

The  solutions  (vj , . . . v^)  and  (v£+1, 
constructed  fro*  functions  of  the  for*  e<x,  where  the  valnes  <  are 
eigenvalues  of  Mj  and  Mj,  respectively.  (See  Section  5.2)  The 
solutions  Vj  of  (5.1) (a)  thus  correspond  to  solutions  of  the  partial 
differential  equation  (3.1)  which  are  made  up  of  modes 


(5.4) 


(<i  +  iw*y  +  st 


If  1  ^  j  i  then  Re  k  <0.  Otherwise,  Re  <  >  0.  The  labelling 
process  of  Proposition  4.2  suggests  that  the  modes  (5.4)  may  be  associated 
with  incoming  waves  if  Re  k  <  0  and  outgoing  waves  if  Re  <  >  0.  Ve  may 
therefore  think  of  V£,...V£  as  "incoming'*  components  of  the  solution 
and  Vj£+j,  ...va  as  "outgoing"  components.  This  identification  is 
developed  in  greater  detail  in  Section  5.3.  The  case  of  tangential 
incidence  is  included  in  this  discussion. 

The  blocks  in  (5.2)  are  not  determined  uniquely,  since  one  can  perform 
similarity  transformations  on  Mj  and  Mj  individually.  These  would 
amount  to  linear  changes  of  variables  among  (v1,...v£)  and  •••▼„) . 

Such  transformations  do  not  alter  the  classes  of  incoming  and  outgoing 
solutions. 

Ve  now  transform  the  boundary  condition  (5.1) (b). 


[I,-S]u(0,u, s)  -  gU.s), 


Since  v  “  Q  *u,  this  can  be  written  as 


(5.5) 


II»~SJQv 


(for 


0) 


'MW'i 
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Let  4i«**«4a  denote  the  column*  of  Q,  end  let 


Q1  -  [*!....«£] 

Q11  “  •  •  •<lnI 

The  boundary  condition  (5.5)  is  thns 

[I.-SHQ1,®11]*  -  g, 
or 

[I»-S]QIvI  -  -[I.-SIQ11?11  +  g 
This  will  be  denoted  by 

(5.6)  N(»,s)v*  *  PCm.sIv*1  +  g  (for  x  *  0) 

It  is  now  possible  to  see  that  the  U.K.C.  can  be  regarded  as  a 
solvability  condition.  First  suppose  that  are  orthonormal, 

so  that  N(u,s)  is  a  matrix  N(», s)  of  the  form  (3.19).  In 
this  case  the  U.K.C.  says  explicitly  that  Idet  N(u,s)l  2.  6  for  Re  a  >  0, 
w  e  IR®.  The  linear  system  (5.6)  can  thus  be  solved  for 
v1,  i.e.,  we  can  solve  for  the  "incoming"  variables  ^....v^  in 
terms  of  the  "outgoing"  variables  . .  .vQ  and  the  boundary 

data  g . 

The  same  conclusion  can  be  reached,  in  at  least  many  cases,  if 

are  not  orthonormal.  Some  arguments  given  in  Section  3.3  show 
that  if  the  qj(u,s)  are  piecewise  continuous  in  (u, s)  and 
homogeneous  of  degree  zero,  then  the  U.K.C.  implies  that  Idet  N(m,s)| 
is  bounded  away  from  zero.  The  linear  system  (5.6)  is  thus  solvable. 

In  the  following  discussions  we  assume  that  the  vectors  qj(«, s)  are 
homogeneous.  This  is  not  a  major  restriction,  since  any  scaling  of  the 
variables  (w,s)  in  (5.2)  can  be  confined  to  M  and  the  blocks  Mj 
and  Mj.  Also  note  comment  (3)  in  Section  3.3. 

The  solvability  condition  contains  a  uniformity  which  is  of  interest. 
Let  s  *  q  +  i£,  fix  q  >  0,  and  consider  the  limit  l«el^  +  1^1^  -+  •. 
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Because  Idet  N(w, s)l  2  $  for  all  (»,  s)  with  Re  s  >  0,  the 

■atria  N(w, s)  is  uniformly  invertible  as  the  frequencies  tend  to  infinity. 

The  significance  of  this  becomes  apparent  when  the  D.E.C.  is  conpared 
with  the  requirement  (3.18),  det  N(u,s)  #  0  for  Re  s  >  0,  which  was 
shown  to  be  necessary  for  well-posedness.  For  reasons  of  homogeneity, 
N(w,s)  »  N(u'.s’),  where 


(Ul2  +  lsl2)1/2 

(See  Figure  5.1)  The  condition  (3.18)  requires  det  N(u',s*)  ^  0  for 
Re  s'  >  0,  but  it  allows  the  possibility  that  det  N(m',s')  may  tend  to 
zero  as  Re  s'  -*■  0. 


Figure  5.1 

That  is,  for  any  fixed  q  >  0,  N(u,s)  may  become  more  and  more  singular 
as  l«o  I2  +  I  Cl2  "*■  ®.  This  may  have  major  effects  on  the 
solution  v2  of  the  system  (5.6):  unless  the  right-hand  side  of  (5.6) 
satisfies  special  constraints,  the  "incoming"  components  v2  may  be 
large  relative  to  v22  and  g  for  large  frequencies.  This  corresponds 
to  a  loss  of  derivatives  at  the  boundary,  or  "weak"  well-posedness.  This 
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will  bo  discussed  ia  greater  detail  in  Section  5.4. 

Ia  Section  3.1  we  assumed  that  the  boundary  condition  has  the  special 
form  (3.4), 

u1  “  SuXI  +  g  (for  x  *  0) 

The  arguments  given  in  the  present  section  make  it  possible  to  justify  this 
assumption.  Suppose  that  we  use  a  more  general  linear  boundary  condtion 

(5.7)  Bu(0,y,t)  ■  h(y,t), 

where  B  is  a  constant  matrix.  Apply  Fourier  and  Laplace  transforms  with 
respect  to  y  and  t,  respectively,  and  write  the  result  as 

BIuI(0,»,s)  +  BXIuXI  »  h(u,s) 

Here  Bx  and  BXX  are  matrices  whose  columns  are  the  first  &  and 
last  n  -  1  columns  of  B,  respectively. 

Now  suppose  w  *  0,  i.e.,  consider  waves  having  phase  velocities 
which  are  normal  to  the  boundary.  In  this  case  N(m,s)  “  *A-X.  Since 
A  is  diagonal  and  has  the  form  (3.2),  one  can  identify  u1  as  the 
"  incoming* '  portion  of  the  solution.  If  we  are  to  be  able  to  solve  for 
u1  in  terms  of  uXI  and  h,  then  B1  must  be  an  invertible 
square  matrix.  Thus  Bx  is  &  x  &,  and  (5.7)  can  be  written  as 

u1  -  -(BI)”*Biru11  +  (BX)_1h. 

This  has  the  form  (3.4). 

5.2.  Solutions  of  the  system  (5.1). 

In  the  previous  sub-section  we  used  modes  (5.4)  for  which  Be  <  >  0, 
i.e.,  which  grow  exponentially  as  x  increases.  This  may  appear  totally 
unreasonable,  since  these  modes  do  not  have  finite  energy  on  the  interval 


0  <  x  <  •.  Ia  addition,  none  of  the  node*  (5.4)  satisfy  the  homogeneous 
initial  condition  which  was  nsed  to  obtain  (4.2)(a)  (  or  (5.2)(a>).  In 
order  to  aiake  the  nse  of  snch  nodes  seen  a  little  nore  legitimate,  we 
therefore  derive  and  discuss  some  representations  of  solutions  of  (5.1).  It 
will  becone  clear  that  the  U.K.C.  has  a  natural  interpretation  in  terns  of  a 
two-point  boundary  value  problem  associated  with  (5.1). 

We  first  justify  a  renark  nade  at  the  beginning  of  Section  4.  There  we 
assumed  that  the  effects  of  the  initial  data  f  in  the  IBVP 
(3.1), (3. 3), (3. 4)  can  be  absorbed  into  the  forcing  tern  F.  This  is  done  in 
order  to  facilitate  the  nse  of  the  Laplace  transform  with  respect  to  t. 

The  prelininary  transformation  can  be  accomplished  by  a  procedure  used  by 
Rauch  [21]  in  a  study  of  the  regularity  of  solutions  of  the  IBVP. 

Suppose  that  the  initial  value  f  vanishes  identically  in  a 
neighborhood  of  the  space-tine  corner  *  *  0,  t  ■  0;  and  suppose  that  the 
pure  initial-value  problem  for  (3.1)  is  well-posed.  There  then  exists  a 
function  w  such  that 

r  2s 
wt  "  Awx  +  L  Bj  j 

w(x.y.O)  *  f(x,y) 

(le  are  assuming  C  ■  0) .  Because  of  the  finite  propagation  speed  associated 
with  hyperbolic  systems,  there  is  a  6  >  0  such  that  w(0.y,t)  =  0  if 
t  j  t,  Choose  h  e  C“ OR)  so  that  h(0)  »  1  and  h(t)  *  0  for  t  >.  6. 

Let  .u  be  a  solution  of  (3.1) , (3.3) , (3.4) ;  and  let  v(x,y,t)  ■  u  -  h(t)w. 

It  follows  that  v  satisfies 

vt  "  Avx  +  l  Bj  +  <F  -  k'(t)w) 

▼  (x,y,0)  »  0 

v*  *  Sv11  +  g,  for  x  =  0. 

This  is  the  form  which  was  used  to  obtain  (4.2).  By  a  re-labelling,  we  may 
continue  to  denote  the  forcing  term  by  F(x,y,t)  and  the  solution  by 
u(x,y,t).  Because  of  the  finite  propagation  speed,  we  may  assume  that 
F(x,y,t)  has  compact  support  with  respect  to  x  and  y. 
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The  system  (5.1) (•)  can  be  solved  by  studying  the  block  version  (5.3). 
Ve  assume  that  Mj  and  M2  are  themselves  block  diagonal, 
where  each  block  is  triangular  and  is  associated  with  a  single  eigenvalue. 
If  necessary,  and  M2  can  be  transformed  to  this  form  by  linear 

changes  of  variables  among  Oj,...V£  and  v$,+i»  •••▼,!•  respectively. 

Under  these  assumptions  each  equation  in  (5.3)  has  the  form 


(5.8)  or 


dw 

—  «  <w  +  H(x) 
dx 

dW  r< 

—  ■  <w  +  L  c.i.  +  H 

dx  J  J 


In  (5.8) (b)  the  Wj  are  solutions  of  equations  which  can  be  solved 
independently  of  (5.8) (b).  We  seek  solutions  of  (5.4)  which  are  in 
L^(0,®)>  since  we  are  ultimately  considering  solutions  of  the  partial 
differential  equation  which  satisfy  energy  estimates  like  (3.5). 

We  first  consider  the  case  Re  <  <  0.  When  one  multiplies  (5.8) (a)  by 
the  integrating  factor  exp(-Kx)  and  integrates,  the  result  is 

x 

(5.9)  w(x)  -  eK*w(0)  +  /  eK(x'z>H(*)dz. 

0 

o 

The  solution  w  is  in  L  (0,«),  since  Re  <  <  0  and  B  has  compact 
support.  The  first  term  in  (5.9)  represents  the  propagation  of  the  initial 
data  w(0);  and  the  second  can  be  regarded  as  a  superposition  of  pulses, 
each  of  which  appears  at  a  point  x  and  is  then  propagated  by  the  natural 
frequency  in  the  problem.  The  more  complicated  equation  (5.8) (b)  can  be 
treated  in  a  manner  which  is  similar  to  the  above. 

When  Re  <  >  0,  the  representation  (5.9)  is  not  appropriate.  In  this 
case  the  function  w  in  (5.9)  would  be  in  L^(0,«)  only  if  the 
initial  value  w(0)  and  the  forcing  term  B  satisfy  a  special  relation. 
Instead,  one  should  impose  the  conditicr  w(x)  ->-0  as  x  •>  «®.  Under  this 
assumption  the  solution  of  (5.8) (a)  has  the  representation 


(5.10) 


/*  eK(x~*)E(z)dz 
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An  analogous  formula  can  be  found  for  (5.8) (b). 

The  system  (5.3)  (or  (5.1) (a))  should  thus  be  associated  with  a  two-point 
boundary  value  problem,  if  the  solution  is  to  be  in  L^(0,®).  The 
components  v*.  corresponding  to  Re  <  <  0,  should  be  prescribed  at  x  =  0; 
and  the  components  v**.  corresponding  to  Re  <  >  0.  should  be  set  to 
zero  at  infinity. 

Te  now  show  that  these  conclusions  make  sense  in  physical  terms.  The 
solutions  w  in  (5. 9). (5. 10)  can  be  factored  into  the  form 

(5.11)  w(x)  -  e**  c(x), 

for  suitable  c(x).  By  inverting  Fourier  and  Laplace  transforms,  one  can 
see  that  the  solutions  of  the  original  partial  differential  equation  are 
thus  made  up  of  functions 

(5.12)  eKx  +  lw‘y  +  #t  c(x> 

An  inspection  of  the  formulas  (5.9). (5.10)  reveals  that  these  functions  are 
superpositions  of  purses  which  are  propagated  by  modes  of  the  form  exp(<x  + 
i«*y  +  st).  According  to  the  labellirg  process  suggested  by  Propostion 
4.2,  the  cases  Re  <  <  0  and  Re  k  >  0  correspond  to  "incoming" 

(rightward  moving)  waves  and  "outgoing"  (leftward  moving)  waves, 
respectively. 

In  the  two-point  boundary  value  problem  the  * 'leftward  moving"  (Re  <  >  0) 

aTT 

components  v  are  set  to  zero  at  infinity.  This  is  physically 
reasonable,  since  leftward  moving  waves  can  arise  only  when  they  are 
stimulated  by  the  inital  data  or  by  forcing  in  tbe  differential  equation. 

The  process  outlined  at  the  beginning  of  this  sub-section  incorporates  both 
of  these  effects  into  the  "forcing"  term  H  in  (5.8).  In  the  case 
Re  <  >  0  the  envelope  c(x)  in  (5.11)  and  (5.12)  has  the  form 

(5.13)  c (x)  -  /*  e~KXH(z)dz 

n 
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(see  (5.10))  This  envelope  can  be  nonzero  only  when  z  is  within  the 
support  of  H. 

"Rightward  moving"  (Re  <  <  0)  waves  way  be  prescribed  at  the  boundary 
z  “  0,  and  in  the  interior  they  may  also  be  stiaulated  by  forcing  or  by  the 
initial  data.  Both  of  these  effects  are  represented  in  the  foranla  (5.9). 

In  order  to  properly  define  the  "incoming"  coaponents  at  z  ■  0,  it  is 
necessary  to  impose  a  suitable  boundary  condition.  The  U.K.C.  describee  a 
class  of  such  conditions;  if  a  boundary  condition  satisfies  this 
criterion,  then  it  has  the  effect  of  ezpressing  the  "incoming"  aodes  in 
teras  of  "outgoing"  modes  and  the  boundary  data  (the  function  g  in 
(5.1)). 

The  remarks  in  this  sub-section  suggest  an  effect  of  the  labelling 
process  of  Proposition  4.2.  The  cases  Re  <  <  0  and  Re  k  >  0  force  one 
to  aake  particular  choices  for  the  boundary  conditions  for  the  systea  of 
ordinary  differential  equations  (5.1) (a).  The  labelling  guarantees  that 
these  choices  are  physically  reasonable. 

5.3  Structure  of  ' *  incoming* *  and  * 'outgoing* *  aodea;  behavior  at 
tangential  Incidence 


In  Proposition  4.2  and  in  Section  5.1  and  5.2,  we  identified  nodes 
ezp(<z  +  iwy  +  st)  as  "incoming"  or  "outgoing"  portions  of  the 
solution  when  Re  <  <  0  and  Re  k  >0,  respectively.  (Here  Re  s  >  0.)  In 
the  present  sub-section  we  ezaaine  more  closely  the  nature  of  these  aodes 
and  the  validity  of  the  above  labelling.  In  particular,  we  discuss  the 
structure  of  the  modes  corresponding  to  tangential  group  velocity  and  the 
effects  of  the  D.E.C.  on  these  modes.  Ve  also  mention  the  strictly 
decaying  modes  which  occur  when  (w,Q  j  T  .  The  present  discussion 
is  given  in  terms  of  the  linearized  shallow  water  equations  and  the 
linearized  two-dimensional  Euler  equations. 

In  this  sub-section  we  consider  only  the  solutions  of  the  hoaogeneous 

_  d  U 

system  (4.7),  ut  «  Aux  +  L  ®j  ay  •  In  the  previous  sub-section  we  noted  that 
these  solutions  can  be  used  to  build  up  solutions  of  the  aore  general 
inhomogeneous  systea  (3.1)  (with  C  ■  0). 


The  idea  behind  Proposition  4.2  is  to  pertnrb  an  oscillatory  solution 
exp  (ierx  +  iw*y  +  i£t)  of  (4.7)  to  a  solution  of  the  foni 

(5.14)  e**  +  iw'y  +  xt 

where  K  *  p  +  io  and  s  =  ij  +  iC*  The  sign  of  Re  K  .  when  Re  s  >  0, 
indicates  whether  the  node  (5.14)  is  "incoming"  or  "outgoing".  The 
Proposition  requires  that  <  and  s  be  invertible  analytic  functions  of 
each  other.  As  noted  earlier,  this  exclndes  the  case  of  tangential 
incidence.  Ve  discuss  such  modes  below,  but  we  first  consider  the  analytic 
case. 

The  mode  (5.14)  can  be  written  as 

(5.15)  epx  +  nt  eiax  +  iw*y  + 

We  regard  this  as  an  oscillatory  mode  which  is  modnlated  by  the  amplitude 
function  A(x.t)  =  exp  (px  +  nt). 


Proposition  5.1.  If  the  oscillatory  part  of  (5.15)  corresponds  to 
group  velocity  which  is  not  tangent  to  the  boundary,  then  (5.15)  approaches 
the  configuration 

-v<x~vt>  eiox  +  iw*y  +  i6t 

(5.16)  e 

as  Ul  +  I  £|  -*•  ®  for  fixed  q  >  0.  Here  v  is  the  x-component 
(i.e.,  normal  component)  of  the  group  velocity. 

Proofs  The  exponent  in  the  factor  A(x,t)  in  (5.15)  can  be  written  as 

B-x  +  t)  fi.  Q.'  £ 

q(q  •  The  ratio  ^  is  equal  to  where  p'  »|^| 

and  n'  *|f|#  1»m*  C  *  (w.C)  end  k  I2  -  l»l2  +  £2.  (This  projection 

is  similar  to  that  of  Figure  5.1.) 

Since  p'  *  0  when  q*  ■  0  in  the  present  case, 


p '  Ap  *  dp  * 

n'  An*  dn' 


n’-o 


The  approximation  becomes  accurate  is  the  limit  n'  0,  i.e.,  when  n  >  0 
Is  fixed  end  lul  +  ICl  ".  Because  of  the  analytic  dependence, 
dp  *  I  dK  do  1 

«*n'  ln*.o  d#  d^  v 

(See  (4.22))  The  exponent  in  A(x,t)  is  thns 


t»(^x  +  t)  *  n("  y*  +  t)  *  -  y(x  -  vt) 


The  Proposition  then  follows. 


The  mode  (5.16)  is  an  oscillatory  wave  which  is  modnlated  by  an 
exponential  envelope  which  moves  at  the  normal  component  of  group  velocity. 
If  v  >  0,  then  (5.16)  decays  as  x  increases.  As  the  group  velocity 
approaches  tangential  incidence,  i.e.,  as  v  •+•  0,  the  mode  decays  more  and 
more  rapidly,  so  that  the  effects  of  the  mode  are  concentrated  near  the 
boundary. 

If  v  <  0  (i.e..  Re  <  >  0),  then  the  modes  (5.14),  (5.15)  are 
"outgoing",  i.e.,  leftward  moving.  According  to  the  remarks  in  Section 
5.2,  modes  of  this  type  are  stimulated  in  the  interior  by  forcing  and  the 
effects  of  initial  data.  The  effects  of  such  stimulation  decay  rapidly  to 
the  left  when  v  <  0  and  Ivl  is  small  (cf.  (5.13)). 

We  now  consider  the  case  where  the  oscillatory  part  of  (5.14),  (5.15) 
has  group  velocity  tangent  to  the  boundary,  i.e,,  v  *■  0. 

The  discussion  in  Section  5.1  suggests  the  effects  of  the  U.K.C.  on 
these  tangential  modes.  In  the  earlier  discussions  the  modes  (5.14)  (5.15) 
have  been  regarded  as  "incoming"  when  p  ■  Re  <  <0.  According  to  the 
remarks  in  Section  5.1,  the  U.K.C.  says  that  the  boundary  condition  must 
give  the  values  of  the  portion  of  the  solution  corresponding  to  Re  K  <  0. 

In  the  present  case,  this  means  that  the  boundary  condition  must  prescribe 
the  behavior  of  waves  moving  tangent  to  the  boundary.  It  is  therefore  of 
interest  to  study  the  structure  of  such  modes.  The  main  conclusion. 
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Proposition  5.2,  will  be  useful  in  the  discussion  of  weak  well-posedness 
given  in  Section  S.4. 

In  the  following  discussion  of  wave  propagation  near  tangential 
incidence,  we  consider  only  the  case  of  constant  coefficients.  When  the 
coefficients  are  variable,  major  complication  can  arise.  For  example,  the 
Hamilton-Jacobi  equations  (4.25)  imply  that  saall  variations  in  the 
coefficients  can  cause  saall  changes  in  frequency  for  a  given  mode.  If  a 
point  («r(v) ,«(r) ,  £(r) )  is  confined  to  ft  ,  and  if  («*,£)  is  near 
the  edge  of  T ,  then  the  aode  can  switch  quickly  froa  "incoming'*  to 
"outgoing"  as  the  paraaenter  t  is  varied. 

Substantial  work  has  been  done  on  the  propagation  of  singularites 
(i.e.,  high  frequencies)  near  tangential  incidence  for  various  types  of 
equations.  The  existing  theory  is  quite  complicated.  See,  e.g.,  Taylor 
[26]  and  the  references  given  therein. 

We  now  discuss  the  form  of  (5.14) , (5.15) .  Suppose  that  (or, <■>.£)  * 
(o0,»0,50)  is  on  ft  and  («0,50)  lie*  on  the  edge  of  r  ,  so  that 
the  oscillatory  part  of  (5 .14) . (5 .15)  corresponds  to  tangential  group  velocity, 
(cf.  Figures  4.1  and  4.2)  The  arguments  of  Propositions  4.2  and  5.1  are 
not  valid  in  a  neighborhood  of  such  a  point,  since  k  and  s  are  not 
analytic  functions  of  each  other.  In  the  case  of  the  shallow  water  and  Euler 
equationa,  there  is  a  square  root  singularity  in  k  as  a  function  of  s. 
(cf.(4.17))  This  is  suggested  by  Figure  5.2;  this  is  a  cross-section  of  the 
cone  ft  corresponding  to  fixed  u  *  wQ.  In  particular,  the  real 
part!  of  <  and  s  satisfy  a  relation 

p  58  R/tj 

for  small  q  >  0.  Here  R  is  a  constant. 
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£  *  Ira  s 


cf  -  Im  k 


Figure  5.2.  Cross-section  of  Q  for  fixed  «  “  «0. 


Proaoeition  5.2.  Suppose  that  the  oscillatory  pert  of  (5.15) 
corresponds  to  tangential  group  velocity  (i.e.,  ▼  ■  0),  and  suppose  that  the 
singularity  between  <  and  a  has  the  fora  described  above.  The  envelope 
A(x,t)  in  (5.15)  is  then  approxiaately  equal  to 


(5.17)  exp  (E  (x  +  t)] 


when  the  ratio  iri  is  snail. 


Proof .  As  in  Proposition  5.1,  the  envelope  A(x,t)  has  the  fora 


exp  [q(^*  x  +  t)] 


In  the  present  case  p*  ~  R/rp  when  q'  is  saall.  Since  q* 
A(x,t)  can  be  written  as 


A(x,t)  ■  exp  *  +  t)J 

*  exp  (Ry<iTcl  x  +  qt] 

This  is  equivalent  to  (5.17).  This  completes  the  Proposition. 

Now  suppose  E  <  0,  i.e.,  Ee  k  <  0.  In  this  case  the  node  (5.15), 


A(x,t)  exp  (iox  +  iet'y  +  lit), 

decay*  rapidly  away  fro*  the  boundary.  Proa  (5.17)  it  is  apparent  that  the 
rate  of  decay  increases  without  bound  as  Id  ■*  m  with  q  >  0  fixed. 

This  situation  contrasts  with  the  earlier  case  v  d  0,  where  the  Halting 
configuration  (5. Id)  is  reached. 

The  behavior  in  the  present  case  corresponds  to  results  froa  the  theory 
of  propagation  of  singularities.  In  the  latter  case  it  is  known  that  as  the 
frequencies  tend  to  infinity,  the  corresponding  portions  of  the  solution 
follow  the  bicharacteristic  curves  (see  (4.25))  better  and  better.  In  the 
present  case  the  directions  of  propagation  lie  along  the  boundary. 

Up  to  now  we  have  considered  nodes  (5.14),  exp  (<x  +  iw*y  +  st)  which 
are  obtained  by  perturbing  purely  oscillatory  nodes  exp(i<rx  +  iwy  +  i£t). 
Such  aodes  correspond  to  (a,£)  s  T  or  points  (o,u,£)  on  the  plane 
in  Figure  4.1.  However,  there  exist  nodes  (5.14)  which  do  not  fit  this 
description.  These  correspond  to  (a,£)  in  the  exterior  of  F  ,  and  for 
Re  s  ■  0  they  have  the  fora 

(5.18)  e<x  +  +  iCt 

One  aay  object  to  labelling  such  aodes  as  "incoaing"  or  "outgoing" 
according  to  whether  Re  <  <  0  or  Re  k  >  0,  since  they  do  not  correspond 
to  the  oscillatory  propagating  waves  aentioned  above.  However,  the  U.K.C. 
still  requires  that  the  boundary  conditions  prescribe  the  behavior  of  such 
aodes  when  Re  <  <0. 

5.4  Weak  RplkrmtflMll 

In  this  sub-section  we  discuss  soae  weak  foras  of  well-posedness  which 
can  be  encountered  when  the  U.K.C.  is  not  quite  satisfied. 

In  Section  5.1  we  considered  the  boundary  condition  (5.6), 

N(», s)vX  -  P(w,s)vn  +  g 
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The  U.I.C.  (3.20)  implies  that  this  system  is  uniformly  solvable  for  v* 
as  Ul  +  I  ^1  -*■  ®  with  i)  >  0  fixed  (cf.  Figure  (5.1)).  At  the  end 
of  Section  5.1  we  remarked  that  the  weaker  necesssary  condition  (3.18), 
det  N(w, s)  ft  0  for  Re  s  >  0,  allows  the  possibility  that  det  N(w,s) 

may  tend  to  zero  as  Re  s  ^  0.  In  cases  where  this  oocnrs,  the  left  side 

of  (5.6)  becomes  singular  as  l<u|  +  |£  |  with  q  >  0  fixed,  so  that 
v1  could  be  large  relative  to  v**  and  g  at  large  frequencies. 

The  incoming  components  wonld  thus  be  less  smooth  than  g  and  the  other 
components,  i.e.,  we  "lose  derivatives"  at  the  boundary.  In  energy 
estimates  like  (3.5),  derivatives  of  data  wonld  appear  on  the  right  side. 

Whether  this  phenomenon  actually  occnrs  depends  on  the  strnctnre  of  the 
right  side  of  (5.6).  For  example,  it  wonld  not  occnr  if  g  **  0  and  if 
P(«, sjv1*  stays  within  the  range  of  N(«, a)  in  some  suitable 
sense.  An  appropriate  characterization  is  that  the  ratios  in  Cramer's  rnle 
stay  bounded  as  Re  1  ^  0.  However,  the  loss  of  derivatives  must  occur  if 

PU.slv11  moves  out  of  the  range  of  N(»,  s)  or  if  arbitrary  g 

are  considered.  The  degree  of  the  derivatives  which  are  lost  depends  on  the 
order  of  the  pole  in  Cramer's  rale.  An  example  will  be  discussed  at  the  end 
of  this  sub-section.  (See  Proposition  5.3) 

The  above  phenomena  have  a  physical  interpretation.  Suppose  g  *  0,  so 
that  the  boundary  condition  (5.6)  is  a  reflection  condition,  i.e.,  the 
"incoming"  modes  v*  are  reflections  of  the  "outgoing"  components 
v11.  The  nature  of  the  reflection  is  governed  by  the  structure  of  the 
boundary  condition.  When  derivatives  are  lost,  the  amplitudes  of  the 
reflected  and  incident  waves  have  ratios  which  tend  to  infinity  as  lu|  + 

Id  -*■  *».  That  is,  the  reflection  coefficients  tend  to  infinity.  In  the 
other  case  mentioned  above  the  reflection  coefficients  remain  finite. 

The  loss  of  derivatives  can  cause  particular  difficnlties  when  the 
spatial  region  is  a  bounded  domain  rather  than  a  half-space.  In  this  case 
waves  can  reflect  back  and  forth  between  various  portions  of  the  boundary, 
so  that  more  and  more  derivatives  are  lost  as  time  progresses. 

Trefethen  [27],  [29]  has  studied  phenomena  analogous  to  the  above  in 
connection  with  the  stability  of  finite  difference  approximations.  He 
emphasizes  the  distinction  between  the  cases  of  finite  and  infinite 


reflection  coefficients,  end  he  describes  how  these  effect  the  precticel 
netore  of  difference  aethods.  In  perticnler,  he  obteins  growth  rete 
estiaetes  for  sm  mild  fores  of  instability. 

The  cese  of  week  well~posedness  is  soaetiaes  essocieted  with  the  tera 
* ' generelized  eigenvalue' 1 .  Recell  the  eigenvelne  problea  (3.12), 


(5.19) 


,  d4  r 

(e)  s4(x)  -  A  dx  +  (1  iwjBj)t  ;  x  >  0 

(b)  4I(0)  -  S4JI(0)  j 

(c)  4  e  L2(0,«») 


If  det  N(w,s)  *  0  for  soae  (w,e)  with  Re  s  ■  0,  then  e  solution 
4  of  (5.19) (e)(b)  exists.  However,  it  aay  feil  to  setisfy  the  boundery 
condition  et  infinity,  4  e  L2(0,*»),  since  the  eigenvelues  of 
M(»,s)  cen  be  purely  iaeginery  when  Re  s  ■  0.  (See  Section  4.3)  If 
this  is  the  cese,  then  the  valne  s  in  (5.19) (e)  is  seid  to  be  e 
’ 'generelized  eigenvelne.''  The  solution  4  does  heppen  to  be  in 
L2(0,-)  if  (e,Q  s  r  end  if  4  is  essocieted  with  the 
eigenvelne  k  of  M(w,s)  for  which  Re  K  <0.  (cf.(5.18)) 

Te  now  deecribe  soae  exeaples  of  week  well-posedness.  We  aeinly 
consider  soae  exeaples  which  were  discussed  by  Kreiss  [11]  froa  a  point  of 
view  which  is  different  froa  the  one  expressed  above.  He  studied  the  system 


for  i  >  0,  y  i  1,  t  >  0,  with  boundery  conditions 

(5.21)  u^O.y.t)  ■  on11  +  g(y,t). 

Here  a  is  e  complex  constant. 

A  ehort  calculation  shows  that  the  eigenvelues  k  of  M(w, s)  setisfy 

(5.22)  k2  ■  s2  +  a2 
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1) 


For  each  dete ruination  of  <  ,  the  corresponding  eigenvector  is  (,  +K, 
Let  k0  denote  the  value  of  <  which  has  positive  real  part*  so  that 
~k0  corresponds  to  the  * ' incoming' *  node.  It  then  follows  that 


(5.23)  N(», s)  - 


(See  the  derivation  of  (5.6).)  If  N(tt,  s)  ”  0.  then 

a2  +  1 

(5.24)  s  -  i  (  2o  )  « 


Kreiss  discusses  three  cases  where  N(t»,  s)  *  0  for  Ee  s  »  0: 

(1)  o  -  +  1 

(2)  I a  I  >  1,  a  real 

(3)  I a I  ■  1,  a  not  real 

He  constructs  solutions  via  Laplace  and  Fourier  transforms  under  the 
assumption  that  the  initial  condition  is  n(x.y,0)  “  0.  For  ease  (1)  he 
obtains  the  estimate 


(5.25) 


i*  f  f 

;  x  c 


e_2l,t  |u(x,y,t)l2  dt  dx  dy 
e"2nt  (|g|2  +  ||^(2  +  ,^12,  dt  dy> 


where  c  is  a  constant.  He  also  obtains  an  interior  estimate  for  u  in 
terms  of  g  alone;  in  this  case  the  spatial  domain  is  given  by  x  2.  6  >  0. 
y  c  ]R.  This  estimate  is  possible  because  of  the  rapid  decay  in  x  of 
the  modes  which  cause  trouble.  In  case  (2)  an  estimate  like  (5.25)  is 
obtained,  but  is  is  not  possible  to  obtain  a  stronger  interior  estimate.  In 
case  (3)  the  situation  is  the  same  as  for  case  (1). 

In  cases  (1)  and  (2)  the  solutions  6  of  (5.19) (a) (b)  are  not  in 
L2(0,«).  Kreiss  refers  to  the  corresponding  values  of  s  as 
"generalised  eigenvalues  of  the  first  kind"  and  "second  kind", 
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respectively.  Case  (3)  corresponds  to  genuine  eigenvalues  s. 

Ve  nov  interpret  the  shove  resnlts  in  tens  of  idees  which  have  been 
developed  in  the  present  paper.  A  short  calculation  shows  that  the 
characteristic  variety  of  (3.20)  is  defined  by  •  o*  + 
i.e.,  its  graph  is  like  the  one  in  Figure  4.1(a).  except  that  the  plane 
is  not  present. 

In  case  (1).  N(w,s)  ■  0  when  i  ■  +  iw,  i.e.,  *  +  iw. 

(See  (5.24)  and  Figure  5.3)  This  neans  that  <  ■  0  and  that  (w,£)  lies 
on  the  edge  of  T ,  so  that  the  nodes  in  question  correspond  to  tangential 
group  velocity.  The  breakdown  of  the  U.K.C.  at  such  points  sieans  that  the 
boundary  condition  does  not  exert  good  control  over  tangential  waves. 
However,  the  effects  of  such  nodes  are  confined  nainly  to  a  neighborhood  of 
the  boundary.  It  is  therefore  reasonable  to  expect  that  a  weak  estinate 
like  (5.25)  would  be  obtained  for  z  2  0  but  that  a  stronger  estinate 
should  be  possible  for  x  2  t  )  0, 

The  structure  of  the  tangential  nodes  was  described  in  Proposition  5.2 
For  the  case  Re  k  <  0  (''incoming'*)  the  nodes  decay  rapidly  away  fron  the 
boundary,  and  the  rate  of  decay  tends  to  infinity  as  I  as  |  +  |£  I  ->■•>. 

Ereiss  observed  this  kind  of  behavior  in  the  present  exanple. 


Figure  5.3.  Locations  of  points  where  N(«,s)  -  0  in  cases (1) , (2) ,  and  (3) 
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We  now  consider  esse  (2),  where  la I  >  1  with  a  reel.  Here 
N(m, s)  ■  0  at  points  (w,s)  **  (n, i£)  for  which  |£|  >  Ini.  (See 
(5.24)  and  Pignre  5.3)  In  this  case  (n,£)  lies  in  the  interior  of  F, 
so  that  the  U.K.C.  breaks  down  for  frequencies  corresponding  to 
non-tangential  group  velocity.  Because  the  nodes  in  question  can  influence 
the  interior,  one  should  not  expect  to  recover  a  stronger  estiaate  by 
restricting  attention  to  the  proper  subdoaain  x  2  8  >0.  The  structure 
of  the  non-tangential  nodes  was  described  in  Proposition  5.1;  as  lee  I  + 

Ul  “  with  i|  >  0  fixed,  these  nodes  do  not  decay  nore  and  aore 
rapidly,  but  instead  approach  the  Uniting  configuration  (5.16). 

In  case  (3)  we  have  la  I  =  1  with  a  not  real.  Here  |£l  <  |e>|, 
so  that  (w, £)  lies  outside  T.  The  corresponding  nodes  are  the 
strictly  decaying  nodes  nentioned  at  the  end  of  Section  5.2.  (See  (5.18)) 

For  these  nodes.  Re  <  t  0  even  when  Re  s  *  0.  Arguaents  sinilar  to 
those  used  in  Propositions  5.1  and  5.2  show  that  when  Re  <  <0,  the  rate  of 
decay  increases  without  bound  as  ln|  +  l£l  "*■  *  with  tj  >  0  fixed. 

We  on it  the  details.  (Also  see  Kreiss  [11].)  Thus,  as  in  case  (1),  a  weak 
estiaate  like  (5.25)  is  obtained  .'or  i  2  0,  and  a  stronger  estiaate  is 
possible  for  x  2  6  >  0. 

We  now  exanine  whether  the  boundary  condition  (5.21)  for  the  systen 
(5.20)  leads  to  bounded  or  unbounded  reflection  coefficients.  For  this 
exanple  the  transformed  boundary  condition  (5.6)  has  the  fora 


After  a  division  and  soae  simplification,  (5.26)  can  be  written  as 
v1  -  R(«,  s)vIX  + 

where 


(5.27) 


R(w,s) 


2 

iw(s  -  <0)  +  on 
i»(s  +  ko)  +  on^’ 
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and  N(w, s)  is  given  in  (5.23).  In  order  to  obtain  (5.27)  we  used  the 
relation  (5.22).  « *  ■  +  m^. 

R(w,s)  is  the  reflection  coefficient  which  relates  the  " incoming' * 
and  ''outgoing"  components  v*  and  v**.  This  coefficient  does 
not  play  a  role  in  the  example  as  discnssed  by  Ereiss  in  [11];  in  his 
formulation  no  outgoing  waves  can  be  present,  since  u(x,y,0)  *  0  and  there 
is  no  forcing  term. 

Proposition  5.3.  Let  («,£)  be  a  point  such  that  N(w.iC)  *  0. 
and  let  s  *  q  +  i£  .  («  and  £  are  fixed.)  In  case  (1).  |S(m.s)I 
remains  bounded  as  Re  s  0.  In  cases  (2)  and  (3).  lR(w, s)l  tends  to 
infinity. 

Proof.  As  Re  s  +  0,  the  denominator  in  R(w,s)  tends  to  zero,  since 
N(t», s)  0 .  The  only  way  that  R(«, s)  can  remain  bounded  is  for  the 
numerator  also  to  tend  to  zero.  Thus  -imiCQ  «  Iwkq  at  the  point 
in  question.  Since  m  A  0  in  cases  (1),(2),(3)  (cf.  (5.25)),  we  have 
“  0.  The  relation  k?  *  s2  +  ((5.22))  implies  s  «  +  i», 

which  corresponds  to  case  (1).  Thus  case  (1)  is  the  only  circumstance 
in  which  lR(u, a) I  could  remain  bounded. 

In  this  case,  u  ■  ♦  1.  For  s  **  q  +  i£  =  q  +  i»,  R(»,s) 
simplifies  to 


R(«.s)  -  - 


- 


R  +  <0 

According  to  remarks  made  prior  to  Proposition  5.2,  Kq  ~  c /q  where 
c  is  a  constant.  Thus  R(m, s)  1  as  Re  a  -+  0.  This  completes  the 

proof. 


The  first  part  of  the  above  proof  illustrates  the  comment  made  earlier 
about  P(w,s)  v11  remaining  within  the  range  of  N(», s) . 

Te  conclude  by  mentioning  an  example  of  weak  well-posedness  studied  by 
Majda  and  Osher  [14,  p.  628] .  In  this  case  the  eystem  of  equations  is 


Maxwell's  equations.  (See  (4.20))  For  a  class  of  boundary  conditions  which 
includes  the  perfect  conductor  boundary  condition,  they  find  that  the  D.E.C. 
fails  at  frequencies  corresponding  to  tangential  incidence.  In  their  energy 
estimate  which  is  analogous  to  (3.5),  they  have  derivatives  of  the  solution 
and  of  g  on  the  left  and  right  sides,  respectively.  The  derivatives  are 
taken  with  respect  to  y  and  t. 


6.  Remarks  on  proof »  of 


We  now  describe  soae  aspects  of  tbe  process  by  which  the  U.K.C.  (3.20) 
is  shown  to  imply  that  the  IBVP  (3 .1) . (3 .3) , (3 .4)  is  well-posed.  This 
discussion  is  not  self-contained;  our  purpose  is  to  point  to  some  features 
of  the  existing  literature  and  relate  them  to  ideas  developed  in  the  present 
paper. 

The  main  step  in  the  proofs  of  well-posedness  is  to  produce  an  &  priori 
"energy  estimate"  like  (3.5).  Such  an  estimate  immediately  yields 
uniqueness  of  solutions  and  continuous  dependence  of  solutions  on  the 
prescribed  data  f,F,  and  g.  Existence  of  solutions  is  shown  via  methods  of 
functional  analysis.  (See,  e.g..  Section  3  of  Majda  and  Osher  [14]) 

The  energy  estimate  is  derived  with  the  aid  of  a  "  symmet riser' ' .  This 
is  a  pseudo-differential  operator  which  has  properties  specified  in,  e.g., 
Kreiss  [10,  p.  291]  and  Majda  and  Osher  [14,  p.639].  Once  the  symmetrizer 
is  constructed  the  energy  estimate  is  obtained  readily,  (e.g.,  [10,  p.  281], 
[14.  p.639]) 

For  problems  having  constant  coefficients,  the  symmetrizer  may  be 
regarded  as  a  smoothly  varying  matrix  function  R(», s)  which  acts  on  the 
transformed  problem  (4.2)  discussed  earlier.  (Here  w  e  JR”,  Re  s  2  0) 

In  the  case  of  variable  coefficients,  one  uses  a  corresponding 
pseudo-differential  operator.  For  simplicity,  we  use  notation  appropriate  for 
the  fdrmer  case. 

We  now  make  some  comments  about  the  construction  of  R(«, s) .  In 
comments  (3)  -  (5)  we  consider  systems  whose  characteristic  variety  has  a 
structure  suggested  by  Figure  4.1.  We  will  denote  by  (<■>', s’)  the  scaled 
variables 

(u,  s) 

(•*',•')  -  (»'.*’  +  i£’>  -  - 

<M2  +  Isl2)1'2 

(See  Figure  5.1). 

(1)  One  first  transforms  M(w,s)  to  simple  block  forms  by  making 
appropriate  changes  of  dependent  variable,  as  in  (5.2),  (5.3).  (See  later 
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consent*.)  In  this  esse  the  transformations  snst  be  smooth  functions  of 
(*>,s),  so  that  the  calculns  of  pseudo-differential  operators  can  be  used. 
One  then  constructs  a  symnetrizer  for  the  simplified  system,  and  at  the  end 
the  effects  of  the  preliminary  transformations  are  incorporated  into  the 
final  symmetrizer.  (Note  equation  (4.S)  in  Najda  and  Osher  [14]  and  the 
equation  after  (4.8)  in  Kreiss  [10])  The  reduction  to  block  forms  receives 
considerable  attention  in  the  theory  (See  Lemma  2.4  in  [10],  Assumption  1.9 
in  [14],  end  Appendix  B  of  Hajda  [13]). 

For  different  (»,  s)  one  say  neel  different  block  forms.  The 
constrnction  is  therefore  done  loca’.ly,  i.e.,  one  constructs  R(», s)  in 
conical  neighborhoods  in  the  («, s)  space  and  then  patches  things  together 
via  a  partition  of  unity.  (See,  e.g.,  equation  (4.6)  in  [14]) 

(2)  The  construction  is  straightforward  when  q'  2.  q0*»  lor 
any  fixed  q0*  >  0.  (See  Section  4  of  [10].)  In  this  case  the 

process  resembles  a  method  for  constructing  Liapunov  functions  for  studying 
the  stability  of  nonlinear  autonomous  systems  of  ordinary  differential 
equations  (e.g.,  John  [8]).  There  thus  remains  only  the  situation  where 
V  is  in  a  neighborhood  of  zero. 

(3)  If  («,£)  is  in  the  interior  of  T  ,  and  q*  is  near  zero, 
then  M(«, s)  can  be  diagonalized  (in  many  cases).  One  can  check  that  the 
construction  is  very  straightforward  in  such  a  situation. 

When  q*  »  0,  the  diagonalizability  is  given  by  Proposition  4.1  of 
the  present  paper.  We  thus  consider  q*  4  0.  For  strictly  hyperbolic 
systems  the  eigenspaces  Ej  in  Proposition  4.1  are  all 
one-dimensional,  so  that  M(w, i£)  has  distinct  eigenvalues.  Locally,  one 
then  has  a  complete  set  of  eigenvectors  which  depend  analytically  on 
parameters.  (See  Lemma  2.4  in  [10])  Some  more  general  eases  are  included 
in  the  "block  structure"  assumption  (Assumption  1.9)  of  Wajda  and  Osher 
[14]. 

(4)  For  (w,£)  outside  T  ,  M(«,i£)  has  eigenvalues  k  for  which 

Re  <  +  0.  (See  Section  4.2  of  the  present  paper.)  These  eigenvalues  can 


be  Isolated  into  blocks  end  handled  as  in  consent  (2)  above.  The 
eigenvalues  of  M(u, s)  which  are  purely  imaginary  when  Re  a  *  0  can  be 
treated  as  in  comment  (3). 

(5)  The  main  difficulty  in  the  construction  of  R(«, s)  occurs  in 
neighborhoods  of  points  where  («,£)  lies  on  the  edge  of  T  and  n  ■  0. 

These  points  correspond  to  tangential  group  velocity.  Ve  noted  in  Section 
4.3  that  M(m, i£)  is  defective  in  this  case. 

As  before.  M(w,s)  is  transformed  smoothly  to  block  form.  The 
defective  block  (or  blocks)  is  handled  in  a  complicated  manner  which  is 
described  in  Section  4  of  Kreiss  [10].  The  construction  is  also  surveyed  in 
Appendix  B  of  Majda  [13].  In  the  case  of  the  shallow  water  and  Enler 
equations  there  is  a  single  defective  block  of  dimension  2x2.  One  may 
wish  to  follow  Kreiss'  construction  for  the  2x2  case. 

(4)  During  the  course  of  the  construction  it  is  necessary  to  relate 
the  symmetrizer  to  the  boundary  conditions.  This  is  the  point  in  the  proof 
where  the  tJ.K.C.  is  nsed;  this  assumption  enables  one  to  solve  for 
"incoming"  (Re  <  <  0)  components  in  terms  of  "outgoing"  (Re  tc  >  0) 
components  and  the  boundary  data  g.  (cf.  Section  4  of  Kreiss  [10].) 
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